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NOTATION 


The.  system  of  notation  proposed  in  SNAKE,  Technical  and  Research 
Bulletin  No.  1-5,  Reference  1,  is  used  in  this  report  wherever  possible. 
The  notation  for  second  and  third  partial  derivatives  is  taken  from 
Reference  2, 

C  Stability  criterion 

I  ,1,1  Momenta  of  inertia  about  x,y,x  axes,  respectively 

K,M,N  Rolling,  pitching,  and  yawing  momenta,  respectively 

N  Typical  static  moment  derivative;  derivative  of  a  moment 

component  with  respect  to  a  velocity  component  dH/dv 

N  Typical  third  partial  derivative;  partial  derivative  of  a 

moment  with  respect  to  a  velocity  component  and  to  an 
angular  velocity  component  yx/dvdr* 

K.  Typical  moment  of  inertia  coefficient;  derivative  of  a 

moment  component  with  respect  to  an  acceleration  component 

UP  Ship  length  between  perpendicular!  (used  ce  a  characteristic 

length  of  body  for  condiment ioualiain*  purpose) 

■  Hess  of  body 

nK|i.  Propeller  revolutions  per  second  At  time  t  and  t.,  reapact- 

c  C1  ively.  1 

p,q,r  Angular  velooltlen  of  roll,  pitch,  and  yew,  reapect Ively 

jM*f  Angular  accelerations  of  roll,  pitch,  sad  yaw,  respectively 

rate  gate  of  deflection  of  rudder  or  other  control  surface 

R.  XonletMce 

c 

T  Propeller  thrust 

t  Thrust  deduction  coefficient 

t.  At  Time  aad  time  interval,  respectively 

♦References  are  listed  on  page  79 

v 


t^  Time  lag  la  control  surface  system 

U  Velocity  of  origin  of  body  axes  relative  to  fluid 

u,v,v  Longitudinal,  transverse,  and  normal  components,  respectively, 

of  the  velocity  of  the  origin  of  body  axes  relative  to  fluid 

Velocity  in  initial  equilibrium  condition:  straight  ahead 
motion  at  constant  speed  with  rudder  amidships 

Au  u“ul 

u,v,w  Longitudinal,  transverse,  and  normal  components,  respectively 

of  the  acceleration  of  the  origin  of  body  axes  relative  to 
fluid 

X, Y,Z  Longitudinal,  lateral,  and  normal  components,  respectively, 

of  hydrodynamic  force  on  body 

Typical  rotary  force  derivative;  derivative  of  a  force  com¬ 
ponent  with  respect  to  an  angular  velocity  component  bY/br 

Yrif  Topical  third  partial  derivative;  partial  derivative  of  a 

‘f  force  with  respect  to~an  angular  velocity  component  and  to 

a  rudder  deflection  ©Y/br  wf 

Y.  Typical  inertia  coefficient;  derivative  of  a  force  component 
with  respect  to  an  angular  acceleration  component  b Y/br 

x,y,*  Body  axes  fixed  in  ship;  x.  y,  and  s  positive  forward,  star¬ 

board,  and  downwards,  respectively.  Origin  of  axes  system 
tot  necessarily  at  center  of  gravity 

Coordinate*  of  center  of  mss  relative  to  body  axes 

?«»*«  Coordinates  relative  to  the  fix  i  earth  axes 

U|  v  O 

xoO'yoQ»*oa  Coordinate*  of  origin  of  body  axes  relative  to  the  fixed 

earth  axes 


%»%  fropeiler  fccxque  m  s  ti m  t  and  tv  respectively 

1 

(5  Angle  of  drift 

5  Angular  displacement  of  a  control  surface,  normally  the 

rudder  angle 

0,9,)^  Angles  of  roll,  pitch,  sad  yaw,  respectively 

p  Mass  density 


^t’^2’^3  Roots  of  stability  equation 

A  prime  (')  applied  after  the  symbol  of  a  quantity  indicates  the  nondimen- 
sional  form  of  the  quantity.  The  nondimensional  expressions  follow  SNAME 
nomenclature.  Reference  1. 


ABSTRACT 


This  report  presents  a  computer  program  for  the  solution  of  a  mathe¬ 
matical  model  representing  the  motion  of  a  surface  ship,  giving  predictions 
of  steering  and  maneuvering  qualities.  The  nonlinear  mathematical  model 
based  on  a  third-order  Taylor  expansion  of  forces  and  moments  in  the  equa¬ 
tions  of  motion  is  reviewed .  The  hydrodynamic  force  and  moment  deriva¬ 
tives  representing  the  input  to  the  program  can  be  obtained  from  present 
captive  model  testing  techniques.  Any  motion  of  a  surface  ship  including 
tight  maneuvers  and  loop  phenomenon  recognized  in  the  spiral  maneuver  for 
a  directionally  unstable  ship  should  be  accurately  predictable.  The  com¬ 
puter  program,  which  gives  predictions  for  the  "Standard  Maneuvers,"  turn¬ 
ing  circles,  zig-zag,  and  spiral  maneuver,  is  described,  and  results  of 
sample  calculations  are  included.  Instructions  for  preparation  of  input 
data  for  the  program,  samples  of  the  computer  results,  and  the  FORTRAN 
listing  of  the  computer  program  are  also  given. 


ADMINISTRATIVE  INFORMATION 

Ihe  mathematical  model  and  associated  computer  technique  presentee, 
by  the  author  should  he  considered  as  a  proposal  and  not  the  current 
standard  for  the  David  Taylor  Model  Basin. 


INTRODUCTION 

A  continuous  growth  in  speed  and  size  of  surface  ships,  an  increasing 
density  of  traffic  on  sea  n  ates,  and  the  development  of  sophisticated  con¬ 
trol  systems  for  steering  and  maneuvering  are  some  of  the  factors  which 
have  stimulated  the  quest  for  precisely  establishing  controllability  quali¬ 
ties  inherent  in  a  surface  ship  design.  As  a  result  the  number  of  ships 
for  which  model  steering  and  maneuvering  trials  are  requested  and  carried 
out  during  full-scale  trials  is  increasing. 


The  time  has  passed  when  a  turning  circle  trial  was  considered  suffi¬ 
cient  for  a  determination;  of  handling  qualities.  Today  it  is  generally 
recognised  that  several  types  of  maneuvers  should  be  known  in  order  to  eval¬ 
uate  the  different  modes  of  performance  of  the  ship  ouch  as  steering,  maneu¬ 
vering,  and  turning.  A  set  of  trials  consisting  of  a  35-deg  turning  circle, 

the  20-20  deg  rig-rag,  and  the  spiral  maneuver  have  been  proposed  for  this 
3 

purpose.  These  maneuvers  are  subsequently  referred  to  as  the  "Standard 
Maneuvers." 

An  adherence  to  these  "Standard  Kaneuvers"  in  both  model  test  and  full- 
scale  trials  should  make  it  feasible  to  establish  criteria  for  steering, 
maneuvering,  and  turning,  and  in  the  future  to  evaluate  precisely  t^eee 
qualities  of  ship  designs.  Another  advantage  of  using  "Standard  Maneuvers" 
as  basis  for  criteria  is  that  the  evaluation  of  ship  performances  can  be 
based  on  a  language  that  is  common  to  operators  as  well  as  to  designers  and 
experimenters. 

Different  testing  techniques  are  in  use  at  model  basins  for  establish¬ 
ing  the  steering  and  maneuvering  qualities  of  a  ship  design.  By  far  the 
most  instructive  are  based  on  free-running  models,  the  performance  of  which 
are  obtained,  for  example,  by  a  direct  execution  of  the  "Standard  Maneuvers" 
in  model  scale.  Despite  obvious  advantages  such  as  direct  modeling  of  maneu¬ 
vers,  the  free-running  model  technique  may  present  difficulties  because  of 
troublesome  scaling  laws,  which  hardly  can  be  taken  into  account  in  this 
technique. 

The  technique  advocated  in  this  report  utilizes  captive  model  testing 
for  the  measurement  of  hydrodynamic  derivatives  with  a  successive 
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prediction  of  the  "Standard  Maneuvers"  obtained  from  a  solution  of  the 
equationo  of  notion  by  means  of  a  digital  computer  or  an  analog  computer 
setup . 

Captive  model  testa  are  performed  by  Earns  of  teat  facilifciea  ouch  as 
the  rotating  arm,  oscillators,  and  the  planar  motion  mechanism.  They  have 
in  the  past  been  adopted  primarily  for  the  measurement  of  the  linear  hydro* 
dynamic  force  and  moment  derivatives  necessary  for  establishing  the  inherent 
directional  stability  of  a  chip  design.  Furthermore,  the  hydrodynamic  force 
and  moment  derivatives  have  been  uced  in  combination  with  the  linearised 
equationo  of  motion  for  analysing  the  turning  ability  of  stable  chipa  in 
the  linear  range.  Hcnevet ,  the  linear  theory  mould  not  in  general,  he  applic¬ 
able  for  predictions  of  the  "Standard  llaacuvorn" ,  aa  it  fails  to  predict 
accurately  the  tight  maneuvers  that  moot  chips  arc  capable  of  performing, 
and  it  cannot  predict  the  maneuvers  of  unstable  ships. 

If  the  loop  phenomenon  (recognised  in  the  spiral  maneuver  for  unstable 
ships)  or  the  characteristics  of  tight  maneuvers  have  to  be  accurately  re¬ 
constructed  analytically,  it  is  ncccooary  to  utilise  equations  of  motion 
expanded  to  include  significant  nonlinear  terms  in  the  Taylor  expansion  of 

forces  and  moments.  Such  a  nonlinear  mathematical  model  has  recently  been 

4 

presented  by  Abkowits. 

5  6 

Chiolctt  and  Strom-Tcjcan *  have  adopted  the  nonlinear  mathematical 
model  and  programmed  the  equations  for  a  digital  cespufcar.  On  the  basis 
of  linear  and  nonlinear  hydrodynamic  derivatives  obtained  by  planar  motion 
mechanics  teats,  they  have  computed  predictions  for  the  "Standard  Maneuvers" 
and  demonstrated  the  accuracy  mlth  uhich  maneuvers  can  be  predicted  in  thin 
fashion. 
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The  captive  Model  t eating  technique  has  an  obvious  disadvantage  in  the 
fact  that  no  direct  display  of  the  ship  maneuvers  is  obtained  from  the  model 
test.  If  such  a  display,  however,  can  be  obtained  accurately  using  computer 


programs  or  analog  setups,  this  disadvantage  is  considered  of  minor  impor¬ 
tance.  Tbe  advantages  in  the  technique  are  numerous;  in  particular,  it 
allows  the  experimenter  to  take  scaling  laws  into  proper  account  and,  in  a 
specific  ship  design,  gives  him  a  direct  insight  into  the  factors  which  can 
be  blamed  for  particular  performance  qualities. 

4 

The  nonlinear  mathematical  model  presented  by  Abkowifcx  is  outlined  in 
the  text  which  follows.  2he  equations  have  been  solved  on  a  digital  computer 
programmed  in  FORTRAN  for  the  IBM  7090  at  TKB.  The  program  gives  a  predic¬ 
tion  of  the  "Standard  Maneuvers"  for  surface  ships  on  the  basis  of  hydro- 
dynamic  force  and  moment  derivatives  obtained  from  captive  model  tcata.  The 
computer  program,  designated  as  Applied  Mathematics  Laboratory  (AMI)  Problem 
XFHCl,  is  outlined  and  data  preparation,  result  sheets,  and  graphs,  etc.  are 
described  in  this  report.  Included  also  are  the  results  of  some  sample  cal¬ 
culations,  which  demonstrate  the  usage  of  the  computer  program  and  its 
ability  to  give  detailed  information  with  respect  to  ship  maneuvers.  The 
saaple  calculations  are  primarily  based  on  hydrodynamic  derivatives  for  the 
MARINER  hull  form  published  in  Reference  5* 

The  appendices  include  instructions  for  the  preparation  of  input  data 
and  the  FORTRAN  listing  of  the  program* 
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MATHEMATICAL  MODEL 


The  derivation  of  a  nonlinear  mathematical  model  representing  the 

4 

steering  and  maneuvering  of  a  surface  ship  is  given  by  Abkowitz.  A 
similar  formulation  has  been  used  as  the  basis  for  the  numerical  computa¬ 
tion  in  the  present  computer  program.  For  the  sake  of  completeness  of 
presentation,  the  development  of  the  Abkowitz  mathematical  model  is  out¬ 
lined  briefly;  a  detailed  discussion  can  be  found  in  Reference  4. 

EQUATIONS  OF  MOTION  FOR  A  SHIP  MOVING  IN  THE  HORIZONTAL  PLANE 

A  general  form  of  the  equations  of  motion  for  a  body,  which  is  allowed 
to  move  in  all  the  six  degrees  of  freedom*  is  obtained  with  the  coordinate 
axis  system  fixed  in  the  body  parallel  with  the  principal  axes  of  inertia, 
but  with  an  arbitrary  origin  not  necessarily  at  the  center  of  gravity. 

•ft 

For  this  case  the  equations  are  (see,  e.g«,  References  1  and  4) 

X  «  m  [  u+qw-rv-x0(q2+r2)+yG(pq-f)+*G<pr4-4)] 

Y  ■  m  [  v+ru-pw-yG<r2{-p2)+*G(qr^)+xG<qp+f)] 

Z  »  »[^pv-qu-sG(p2+q2)+xG<rp-4)+yG(rq+^)] 

K  «  Ix$+  (IE»Xy) qr+m  [yG (topv-qu)  -®G  <*fru-pw)j  ^ 

M  «  Iy4+(IX-Ir)rp+M  [  *G  (u+qw-rv)  -x„  (w+pv-qu)] 

H  «  Ixf+(Iy-Ix)pq+»  [  *G  (v+ru-pv)  -y,  (u+qw-rv)] 

wicre  the  left-hand  side  represents  the  forces  and  momenta  along  the  coor¬ 
dinate  axes  and  the  right  hand  side  shows  the  corresponding  dynamic  re¬ 
sponse  terms. 

*lhe  equations  are  developed  assuming  the  mass  of  the  body  as  being  constant 
in  time,  which  can  be  considered  true  for  most  ships. 

5 


When  dealing  with  steering  and  maneuvering  of  our face  chips,  the  primary 
motions  can  be  considered  to  take  place  in  the  horizontal  plane,  and  vertical 
motions  can  ha  noglcgtcd.  Further,  choosing  an  axis  cyotcm  in  the  plana  of 
symmetry  of  the  body  and  assuming  that  the  center  of  gravity  lice  in  the 
centerline  plana  and,  therefore,  y^=0 ,  the  equations  of  motion  for  a  chip 
■roving  in  the  horizontal  plane  hecoae 

X  =  m  [u-rv-xGr  +*npr]| 

Y  =  *  [vf ru-*Gp+xGr] 

K  =  Ixp-a-*G  (vMru) 

K  -  Isf+m*xG(v+ru) 

In  the  following  treatment ,  rolling  and  heel  of  t  ie  chip  has  been 
neglected,  since  thoy  are  felt  to  hava  little  influence  on  steering  and 
maneuvering,  with  the  posoiblo  exception  of  fact  warships,  which  heel 
appreciably  in  turno.  The  cquationo  for  eteoring  and  maneuvering  of  a 
•urfaca  ahip  thus  reduce  to 

X  »  *  [i-tv-x^r2] 

Y  »  «  [skmxgf] 

X  «■  I^f+wXQ^ru) 

TAYLOR  EXPANSION  OP  FORCES  AND  MOMENTS 

The  forces  and  moments  on  the  left-hand  side  of  the  equations  of 
notion  can  bo  enpreoeed  as  functions  of  properties  of  the  body,  properties 
of  the  fluid,  and  properties  of  the  motion.  When  considering  a  specific 
hull  form  and  using  the  generally  accepted  scaling  laws,  the  forces  and 
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moments  ray  be  concidered  as  functionc  of  the  motion  and  orientation 
parameters  only.  When  dealing  x;ith  otcering  and  maneuvering,  they  are 
also  conoid crcd  as  functions  of  the  deflection  6  of  control  curfaceo 
(rudder) : 

Moment  J  ~  f (ProPertia8  of  ootion9  rudder  deflection) 


=  f(xo»yo»*o»^,0,yrsu»v,v,p,q,r,u,v»wsp,4*r,  6  ,6  ,5,  etc.) 

- - - »  “< . — — gB»  «s£ — >■-■» - — ■**» 

orientation  motion  parameters  control 

parameters  surface 

parameters 

When  conoideristg  cation  in  an  unrestricted  horizontal  plane,  it  ic 
clear  that  no  forceo  or  moments  are  exerted  on  the  chip  due  to  a  change  in 
orientation,  and  the  forces  and  cement  will  then  only  be  functions  of  the 
three  degrees  of  freedom  eotion  parameters  and  the  rudder  deflection: 


X  1  i  *i 

Y  f  »  f(u,v,r,u,v,f,6,5 ,0,  etc.) 

N  J 

In  tho  following  troatcont,  it  is  furthor  assumed  that  the  control 

forces  and  cocanto  produced  by  a  dof lection  of  tho  control  curface  (rudder) 

are  duo  to  the  dof  lection  $  only,  while  forces  and  cements  produced  on  the 

•  •« 

ship  as  a  result  of  6  and  6  are  negligible.* 

The  functions  describing  tho  forces  and  cement o  can  bo  developed  into 
a  useful  fora  for  analysis  purposes  by  the  use  of  the  Taylor  expansion  of 
a  function  of  several  variables.  Tho  forces  and  cements  can  thus  be  oh- 
prcoocd  to  any  desired  degree  of  accuracy  by  considering  sufficient  terms 


•  ** 

*Thc  variables  6  and  6  arc  considered  negligible  in  the  treatment  of  ship 
motions,  but  they  are  not  necessarily  negligible,  if  determining  the  forces 
on  the  rudder  Itself j  c.g.,  the  torque  cn  tho  rudder  stock  during  a  maneuver. 


In  the  expansion.  If  the  expansion  ie  lioitcd  to  the  first  order  tern*, 
the  well-known  linearised  expansion  will  bo  obtained. 

If  straight  ahead  cotton  at  constant  speed  with  rudder  amidships  ie 
chosen  as  the  initial  equilibrium  condition,  the  linearised  expansion  of 
the  forces  and  oooent  (Equation  (4)  )  become!* 

X  a  X*  +  XuAu  +  Xyv  +  Xrr  +  X.u  +  X.Y  +  X.f  +  X^6  (5) 

where  An  =  (u-Uj)»  with  similar  expressions  for  Y  and  N. 

Similarly,  the  Taylor  expansion, including  terns  up  to  third  order, 
become® 

x  =  X*  +  [  XuAtt  Xvv  +  Xrr  +  X.u  +  +  X.f  +  X3  &  ] 

+  2  t  W“2  +  W2  +  .  *6«  *2  + 

n-xjn-v  +  2-V»‘r  +  .  *  **  ]  „ 

W 

*  3|  [ +  Vy3  *  •••*•  +  + 

3,5W^  +3'*tauP^v  +  *****  +  3*xr$S*^2+ 

6,Xttv/***Vt  +  6*Xuv^iu*ui)  +  .....  +  6 *X. . ^  vr&  J 

with  alallar  expressions  for  X  and  N* 

LINEAR  MA1HBM4TICAL  mm.  FOR  STEERING  AND  tttNEUVERXKS 

Equating  the  linearised  expansion,  Equation  (5),  with  the  dynamic 
rosponso  terns  given  on  the  right-hand  side  of  the  equations  of  motion, 
Equations  (3),  and  neglecting  dynamic  response  of  cecond-ordcr  smallness 
in  the  same  way  as  second-order  tores  have  been  neglected  in  the  force  and 
moment  expansions,  the  linearized  equations  of  motion  for  steering  and 
maneuvering  are  obtained 
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1 '  S- r  t/PSSSHp1*  - 


X,  4  X  +  X  V  +  X  r  4  X.«  +  X.v  4  X.r  +  Xi  6  = 

*  u  v  r  u  v  r  0 

+  Yu£.u  +  Yyv  4  Yfr  4  Yaa  4  Y^v  4  Y^r  4  Y$6  =  n(vfruj43^,r)  (7) 

U,  4  B  &u  4  N  v  4  Nr  4  N.u  4  N.v  4  II. r  4  Ns  6  =*  I_i+nx_(v+ru.) 

*u  vruvro  Sul 

The  derivatives  Xy,  H.,  Xr»  X.»  and  Xg  arc  ail  aero  for  any  chip  or 

body  with  symmetrical  chape  pare  end  starboard.*  Ac  a  consequence,  Y  ,  Y^, 

4 

N  and  N.  eust  also  be  zero. 
i‘  u 

With  the  teres  ea  the  right-homd  side  of  the  equations  brought  over  to 
the  left  side  end  combined  with  similar  teres,  the  linear  mathematical  model 
for  the  steering  cad  maneuvering  of  a  surface  chip  finally  becomes 

(X.-a)  »u  4  X  £Lu  =0 

u  u 

(\\-a)  -v  4  Yyv  4  (V.-aotG)*f  4  (Yr*» Uj)*r  4  Y^6  «  o  (8) 

(H^-oxg)v  4  Nyv  4  (N^-Xg)  *r  4  (N^aXgU^r  4  Hjfl  «  0 

On  the  basis  of  the  linear  model,  Equations  (8),  the  well -known  criterion 
for  dynamic  otsbility  in  straight  lino  motion  can  be  evaluated  aa 

C  a  Yv(Hr-a«0u1)  -  Hy<Yr-wij)  >  0  (9) 

For  a  dynamically  afcablo  Chip,  the  model  can  furthermore  be  applied 
to  predict  maneuvers  ao  long  as  only  email  rudder  deflections  and  email 
deviationo  from  tha  original  straight  lino  motion  are  considered.  The 
3  imitations  of  the  model  arc,  however,  obvious  from  the  fact  that  no  speed 
loss  is  indicated. 


*Thio  is  one  of  the  advantages  by  choosing  axis  systems  in  the  plane  of 
symmetry  of  the  body. 
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KONLIaEAR  MnTi3H3i\XlCAL  hODEL 

to  obtain  realistic  predictions  of  maneuvers  such  sc  tight  turns  for 
large  rudder  angles  and  to  predict  the  performance  of  a  dynamically  unstable 
ship,  it  becomes  necessary  to  develop  and  eolve  a  nonlinear  mathematical 
model,  which  includes  higher  order  terns  in  the  Taylor  expansion  of  forces 
and  moments. 

The  ncnlincar  mathematical  model  used  as  a  hacia  for  the  computer  prog 
has  been  based  on  a  Taylor  expansion  of  forces  and  moments  including  terms  of 
up  to  third  order j  -sec  Equations  (6) .  The  inclusion  tf  terms  higher  than 
third  order  was  not  considered  to  increase  the  accuracy  of  prediction  sig¬ 
nificantly.  Furthermore,  practical  limitations  of  measurement  techniques 
and  the  state  of  refinement  of  present  theory  did  not  justify  the  inclusion 
of  higher  terms. 

Symmetry  considerations  demonstrate  that  the  X-cquation  should  be  an 
even  function  of  the  parameters  v,  r,  6,  v,  and  r;  similarly,  the  Y-  and 
N-equations  are  odd  functions  of  the  eama  parameters.  Consequently,  odd 
terms  in  v,  r,  6,  </,  and  r  have  been  eliminated  from  the  X-cquation,  and 
oven  terms  in  the  same  parameters  from  the  Y-  and  N-cquationo.  An  alterna¬ 
tive  solution  would  hava  been  to  introduce  absolute  values  of  the  parameters 
v,  r,  6,  v,  and  i  into  the  equations,  but  this  was  considered  less  attrac¬ 
tive. 

As  a  further  consequence  of  the  body  symmetry,  Yu>  Yuu,  Yuuu>  Y. 

and  corresponding  derivatives  in  the  moment  equation  H  ,  N  ,  N  ,  N. 

M  U*  uu*  uuu*  u 

arc  all  zero. 
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An  uncymmetrical  force  (for  inotancc,  the  cide  force  from  a  single 
propeller)  has  been  taken  into  account  by  conotant  teraa  Y*  and  11^.  in  the 
Taylor  expansion.  An  uncymmetrical  aide  force  has  been  considered  a  function 
of  speed,  and  terms  Y^,  H^u,  N^.  have  consequently  been  introduced 

into  the  mathematical  cod  cl  to  facilitate  that  changcG  of  cide  force  with 
speed  are  taken  into  account.* 

The  nonlinear  equationo  can  be  reduced  further  by  conoidering  the 
nature  of  the  acceleration  forceo.  Abkcwite  otatcc,^  that  no  cecond  or 
higher  order  acceleration  tome  can  bo  expected.  This  io  baccd  on  the  accuisp- 
tion  that  there  ia  no  significant  interaction  between  viccouc  and  inertia 
properties  of  the  fluid  and  that  acceleration  forceo  calculated  from  poten¬ 
tial  theory  give  only  linear  terns  when  applied  to  submerged  bodice. 

Abkcwita  further  roaconc  that  terns  repre coating  croon-coupling  be¬ 
tween  acceleration  and  velocity  parcrotorc  arc  aero  or  negligibly  email 
for  reasons  similar  to  those  just  given. 

The  validity  of  these  baoic  considerations  of  Abkcwitc*  hao  been 
vorifiod  by  tho  onporicentai  measurements  scperfced  in  Reference  5. 

Equating  the  nonlinear  Taylor  enpanaion,  Equationo  (6),  with  dynaaic 
raoponco  fcorao,  Equationo  (3),  end  taking  tho  abovo  conoiderafciona  into 
account,  the  nonlinear  equations  of  motion  finally  become 

X-Equation:  (w-X^) u  «  f^(u,v,r,5) 

Y-Equation:  (m-Y^)*  +  (mx^Y^f  =  f2(u,v,r,S)  (10) 

‘•If  an  unoycmotrical  force  chould  turn  cut  to  be  a  function  of  other 
parcmetcrc  than  epeed,  thio  uncymofcry  could  easily  be  introduced  into 
the  prccent  mathematical  model.  It  would  have  been  more  difficult  to 
do  this  if  absolute  values  of  the  parameters  has  been  applied. 
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H-Equation:  (mx_-N.)v  +  (I  N.)r  -  f.(u,v,r,6)  (10) 

G  v  **  r  3  coat'd 


where 


fjCu.v.r.i)  =  X*  ♦  X  Au  +  ixuau2  +  |xuute3  ♦ 

K/+  <KC+”-G)r2+  ¥t&?*  L’*4*  +  KnA1  +  + 

(Xvr+m)vr  -r  Zvp6  +  Xr6r6  +  X^vrAu  +  X^vSAu  +  X^^Au 

f2(u,v,r,5)  =  Y*  +  Y*uAu  +  Y*uyAu2  + 

Y  v  +  rY  v3+  ~Y  vr2+  ~Y  ft v<$2+  Y  vAu  +  ~Y  vAu2+ 

v  6  vw  2  vrr  2  voo  vu  2  vuu 


<Yr-nu) 


r  +  l^rrr^  ^-^Z+  Y  -r^u  +  ^--*,Au‘*+ 


‘rSS1 


TO 


i  j. 

ruu 


V  +  ha5i3+  hrrSr2+  Va“  +  I,rSvli 

f3(u,v,r,6)  =  K*  +  N*uAu+  K*uuAu2  + 

V  +  &A  KrrvA  Ks«v42+  "vuvA"  +  K» A+ 

<W>r  *  £«„/+  ^rA  KjjrA  K^rAu  +  ^saA 

»a4  +  bsu^*  Kv»4y2+  K«4r2+  Vto" +  ^5uu4a“2+  Kvravr4 

PRINCIPLES  FOR  SOLUTION  OF  MATHEMATICAL  MODEL 
USING  DIGITAL  COMPUTER 

METHOD  OF  NUMERICAL  SOLUTION 

The  mathematical  model,  Equations  (10),  can  be  solved  with  respect  to 
the  accelerations  0,  •>,  and  r,  which  become 


fx(u,v,r,6) 


U  a 


V  a 


(»-Xd) 


<VN.)  f2(u,v,rs6)  -  (atXg-Y^)  f3(u,v,r,$) 


(m-Y^) (IS~N. )  -  (m*G-N^) (m*G-Y.) 


(11) 
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(11) 

cont’d 


(“-¥$>  £3(u,v,r,6  )  -  <mxG-iy  £2(u,v,r,4) 

*  =  (m-^)<IE-N.>  -  (n«G-N.) 

These  solutions  can  ba  rewritten  in  the  tom 
gx  [  t,u(t),v(t),r(t),6(t)] 

|f  =  82[  t,\i(t),v(t>,r(t),6(t)]  (12) 

=  S3  [  t  ,u(t)  ,v(t)  ,r(t)  ,6(t)] 

It  is  Geen  tr»s  t  the  mathematical  model  has  been  reduced  to  a  set  of 
three  first-order  differential  equations.  An  approximate  numerical  solu¬ 
tion  for  this  type  of  equations  is  readily  obtained  on  a  digital  computer. 
The  process  in  the  numerical  solution  is  that  the  values  of  u,  v,  and  r  at 
time  t+At  are  obtained  from  knowledge  of  the  values  of  u,  v,  r,  and  &  at 
time  t. 

A  simple  first-order  method  has  been  applied  in  the  computer 
the  valucc  at  time  t+At  are  obtained  simply  by  the  first-order  Taylor 
scries  expansion 

u(t+At)  ■  u(t)  +  At-u(t) 

v(t+At)  »  v(t)  +  At*v(t)  (13) 

r(t+At)  m  r(t)  +  At»r(t) 

This  method  is  found  to  give  adequate  accuracy  for  the  present  type  of 
differential  equations,  because  of  the  fact  that  the  accelerations  u,  v, 
and  f  vary  only  slowly  with  time.  This  is  due  to  the  large  mass  and  inertia 
of  a  ship  compared  to  the  relatively  small  forces  and  moments  produced  by 
its  controTsurfaceaV  < 
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Furthermore,  digital  computers  enable  long  repetitive  calculations 
to  be  made  fast  and  accurately,  and  any  desired  accuracy  of  the  solutions 
can  be  obtained  using  small  time  intervals  At. 


CALCULATION  PROCEDURE  FOR  PREDICTION  OF  TRAJECTORY 


So  far,  the  mathematical  model  has  been  developed  in  dimensional 
form.  The  development  has  on  the  other  hand  been  completely  general,  and 
the  equations  are  equally  valid  in  the  nondimensional  form.* 

In  the  computer  program,  the  mathematical  model  has  been  adopted  in 
its  nondimen aional  form.  To  describe  the  calculation  of  a  trajectory  in 
dimensional  form  on  the  basis  of  the  nondimensional  equations,  the  non- 
dimensional  ized  form  of  a  given  quantity  will  be  indicated  by  the  prime 
of  that  quantity  in  the  following  discussion. 

Assuming  that  a  full  set  of  nondimensional  hydrodynamic  coefficients 
(Xj,  XuuV  YvJ  etc.)  is  available  and  that  the  rudder  deflection  h  is 
defined  as  a  function  of  time,  the  first  step  in  the  calculation  of  the 
trajectory  of  a  ship  maneuver  would  be  to  define  the  initial  condition, 
l.e.,  set  the  nondimensional  values 


u(t)'=  u(t)/u(t) 
v(t)'=  v(t)/u(t) 

(14) 

r(t)'=  r(t)/(u(t)/LBP) 

$(t  )'=  6(t) 

at  time  t=0.  Having  done  this,  the  nondimensional  accelerations  u', 
and  r*  can  be  calculated  from  equations  (11),  and  the  corresponding 
accelerations  in  dimensional  form  from 


*The  velocity  used  for  nondiraensionalization  should  be  the  velocity  at 
any  time,  t  rather  than  the  initial  velocity. 
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(35) 


u(t)  =  u(t)'*  (u(t)2/LBP) 

v(t)  =  v(t)'.  (u(t)2/LBP) 

f(t)  =  r(t)‘.  (u(t)2/LBP2) 


The  new  velocities  in  dimensional  form  at  time  t=At  can  be  obtained  from 
Equations  (13)  and  the  corresponding  nondimensional  values  from  Equations 
(14).  The  process  is  then  repeated  using  the  new  values  for  u',  v',  r', 
and  J*  in  Equations  (11),  and  so  on.  The  values  of  the  velocities  at  a 
time  t  are  thus  obtained  from 


u(t)  =  u(0)  + 

v(t)  =  v(0)  + 


t-At 

u(T)  *At 


t-At 

At 


r(t)  =  r(0)  + 


t-4t 

r(V)  -At 


(16) 


where  u(0),  v(0),  and  r(0)  are  the  values  of  u,  v,  and  r  at  t=0,  and  r 
represents  intermediate  values  of  time  (between  time,  0  and  time,  t-At) 
at  which  the  accelerations  u(t),  v(t),  and  r(t)  are  determined. 

The  instantaneous  coordinates  of  the  path  of  the  origin  of  the  ship 
xoQ(t)  and  yoQ(t)  relative  to  the  fixed  earth  axes,  the  instantaneous 
radius  of  curvature  R(t),  angle  of  yaw  jKt) ,  etc.,  can  be  obtained 
similarly  from  the  velocities  by  using  the  formulas 


t-At 

^(t)  -  f(Q)  +  >  r(Y)  *At 


t-At 


xo0(t):=  xo0(0)+  ^  [  V(T)  *  cos)^+  (u (r) -U(0) )  •  sinyw  )-At 


1*? 


(17) 


t-Afc 

yoO(t>  “  yoO^+  y  [(u (T)-u(O))  •  coa)fiCt)  -  v(T)  *siny(7)  ]  *At 

.  .  _  V; u(t)+u(0))2+  v(t)2  (17) 

*'c'  r(t)  cont'd 

The  accuracy  of  the  predicted  trajectory  can  be  controlled  by 
running  the  calculation  with  different  values  of  the  time  interval  At. 

It  is  found  that  a  high  accuracy  is  easily  obtainable,  and  a  time  inter¬ 
val  of  At  =  1  sec  has  been  chosen  as  standard  in  the  computer  program. 


DEFINITION  OF  RUDDER  DEFLECTION 

It  is  necessary  in  the  calculation  of  a  ship  trajectory , as  mentioned 
above,  to  define  the  rudder  deflection  as  a  function  of  time.  This  has 
been  accomplished  in  the  computer  program  by  assuming  the  rudder  to  move 
with  a  certain  constant  rate  of  deflection  and  assuming  a  certain  timelag 
between  the  instant  the  rudder  deflection  is  ordered,  and  the  instant  the 
rudder  begins  to  mova.  A  rudder  deflection  up  to  a  certain  given  angle 
^const  wou^  ke  executed  in  the  program  as  indicated  in  the  following 
example : 


S(t)  =  6(tx)  until  t  >  tlag+  tx 

then  6(t)  *  +  rate* (t-t^-t^flg)  until  6(t)  *  ^const 

then  j(t)  =  4COMt. 

A  rudder  function  of  this  type  gives  a  close  approximation  to  the 
actual  time  history  of  a  ship’s  rudder  when  a  certain  maneuver  is  ordered 
on  the  bridge,  and  almost  any  practical  rudder  sequence  encountered  when 
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considering  ship  maneuvers  can  be  built  up.  The  zig-zag  maneuver  can,  for 


example,  be  built  up  as  follows,  using  these  principles: 


6(t)  =  6<tL) 

then  6(t)  =  6(tj)  +  rate* (t-t^flg-tp 
— th™  <»>  -  ^const 

then  a<t)  =  iconBt 

than  4(t)  =  acoM£-  rate-Ct-tlag-t2) 

n  then  6(f)  =-6 

■  const 

then  6(t)  =-6  , 

*  “const 

then  6(t)  .+  rate. (t-t,  -t_) 

vv  '  wconst  v  lag  y 

— then  repeat. 


until  t  >  tlag+  tx 

until  S(t)=6  . 

const 

until  t  =  t„  when  ^ff~ b 


const 


until  t  >  tlag+  t2 
until  6{t)5=  -6C 

*3 


const 

until  t  -  when ^=-6 

until  t  >  t 


const 


•i  t0 

lag  3 


until  6(t)=  6 


const 


COEFFICIENTS  IN  MATHEMATICAL  MODEL 


EXPERIMENTAL  TECHNIQUES  FOR  MEASUREMENT  OF  COEFFICIENTS 

To  perform  the  computations  of  ship  maneuvers,  it  is  necessary  to  know 
the  various  hydrodynamic  derivatives  Yy,  N  ,  etc.)  which  appear  in 
the  mathematical  model.  Equations  (10) ,  These  coefficients  depend  largely 
upon  the  ship  geometry  and  design,  and  in  general  they  differ  significantly 
from  one  hull  form  to  another.  For  most  of  the  coefficients,  it  is  necessary 
to  rely  on  model  testing  techniques  of  special  nature  in  order  to  determine 
the  values  for  the  particular  ship  form. 

The  coefficients  are  by  definition  partial  derivatives  of  a  force  or 
moment  with  respect  to  one  or  more  of  the  motion  parameters.  To  obtain  the 
different  coefficients,  it  is  necessary  to  let  the  model  execute  various 
forced  motions  and  to  measure  the  forces  and  moments  as  functions  of  the 
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different  notion  parameters.  An  example  might  illustrate  this  principle. 

For  a  model  which  hao  been  towed  at  different  specific  drift  angles,  corres¬ 
ponding  forces  Y  and  moments  N  have  been  measured.  Figure  1  shows  the  non- 
dimensional  values  Y1  and  N*  plotted  ld  a  function  of  the  nondimensional  side 
velocity  v’=v/u.  From  these  measurements,  it  is  now  possible  to  obtain  the 
derivatives  with  respect  to  the  side  velocity  v,  namely,  Yy,  ^vvv  as  well  as 

N  and  N  .  The  derivatives  are  related  in  a  simple  manner  to  the  coeffi- 
v  vvv  r 

cients  in  the  third-order  poiynoainals,  which  give  the  best  curve  fitting  to 
the  experimental  values.  Thus,  if  the  third-order  polynomials  fitted,  e.g., 
by  a  least  squares  procedure,  are  of  the  form: 


Y*  =  +  a^-v'  +  a^'V*^ 
N*  =  +  b^*v*  +  b^v'^ 


(18) 


then  the  derivatives  would  be  directly  related  to  the  polynom-coeff icients 
as  follows: 


N  *  =  b- 
v  1 


(19) 


Different  testing  facilities  such  as  rotating  arm,  oscillators,  and 
planar  motion  mechanism  are  capable  of  executing  model  tests  with  various 
types  of  forced  motions.  The  most  versatile  instrumentation  is  probably 
the  planar  motion  mechanism  because  any  type  of  motion  with  respect  to  which 
derivatives  are  desired  con  be  produced  by  this  instrumentation.  A  detailed 
discussion  of  a  planar  motion  mechanism  and  the  technique  for  measuring  the 
different  derivatives  for  a  surface  ship  is  presented  in  Reference  6.  Here 
it  is  sufficient  to  mention  that  measuring  techniques  are  available,  which 
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in  model  scale  permit  measuring  the  different  derivatives  appearing  in  the 
mathematical  model.  Equations  (10) . 


CALCULATION  OF  COEFFICIENTS  IN  X-EQUATION 

Three  of  the  coefficients  in  the  X-equation,  X  ,  tX  ,  and  ~X  .  are 

*  u*  2  uu*  6  uuu* 

calculated  in  the  computer  program  on  the  basis  of  the  results  from  open- 
water  propeller  test  and  the  ship  effective  horsepower  data. 

When  the  ship  is  sailing  straight  ahead  with  constant  velocity  u^> 
the  propeller  thrust  working  with  the  thrust  deduction  exactly  equals  the 
resistance  of  the  ship 


X  =  T(l-t)  -  Rfc  -  0 


(20) 


This  equilibrium  condition  defines  the  initial  propeller  thrust  and 
the  corresponding  propeller  torque  and  revolutions. 

As  soon  as  a  maneuver  is  initiated,  this  equilibrium  condition  is  dis¬ 
turbed.  The  X- force,  which  represents  the  difference  between  the  propeller 
thrust  and  the  ship  resistance,  will  vary  as  a  function  of  the  speed. 
Approximating  the  X-force  by  a  third-order  polynomial, 

2  3 

X(u)  s  aQ  +  a^*Au  +  a2  «Au  +  a3  *Au 

where  Au  =  (u-u,),  the  derivatives  X*,  X  ,  ■sX  ,  rX  can  be  obtained 

JL  U  m  UU  O  UUU 

directly  from  the  coefficients  of  the  polynomials  as  follows: 


X*  «  a.c^O;  X  =  a.:  4x  =  a„ :  -rX  -•  a. 

0  *  u  1*  2  uu  2*  6  uuu  3 


1„  .  Ix  _ 

6  uuu 

In  the  program  the  actual  X-force  is  computed  at  the  different  speed 


values  for  which  the  ship  resistance  is  known  from  the  ship  effective  horse¬ 
power  data.  The  corresponding  propeller  thrust  values  are  computed  using 


different  assumptions  that  depend  upon  the  type  of  engine  and  the  engine  net¬ 
ting  to  be  maintained  during  the  maneuver. 

The  propeller  thrust  can  thus  be  calculated,  either  assuming  constant 
propeller  revolutions  or  assuming  the  propeller  torque  to  vary  proportionally 
to  the  revolutions  in  a  certain  power.  If  torque  is  assumed  to  vary  inverse¬ 
ly  proportional  to  propeller  revolutions,  the  thrust  values  corresponding  to 
a  turbine  power  plant  capable  of  maintaining  a  constant  power  output  would 
be  obtained.  If  torque  is  assumed  to  be  constant  during  the  maneuver,  the 
corresponding  condition  for  a  Diesel  power  plant  would  be  obtained. 

SCALE  EFFECTS 

Most  of  the  coefficients  to  be  used  in  the  mathematical  model  would 
be  obtained  from  model  tests,  and  in  this  connection  it  is  reasonable  to 
give  some  considerations  to  scale  effects  in  the  measurement  of  the 
coefficients. 

The  model  tests  would  be  conducted  according  to  Froude's  law,  hence 
the  Reynolds  number  Would  not  be  satisfied,  and  the  possibility  of  Reynolds 
number  effects  should  be  Eecogniied . 

Tests  with  airfoils  covering  a  wide  range  of  Reynolds  numbers  indicate 
that  change  of  Reynolds  number  apparently  has  no  systematic  effect  on  the 
lift-curve  slope.  However,  the  variation  of  maximum  lift  might  be  appreci¬ 
able  because  separation  or  flow  breakdown  occur  earlier  for  the  relatively 
thicker  boundary  layer  around  a  model  body  at  the  lower  Reynolds  number . 

These  results  from  airfoil  testing  can  be  applied  in  the  present  discussion 
of  scale  effects,  as  most  of  the  Y- forces  and  N-moments  would  be  due  to 
similar  lift  and  circulation  effects.  Thus,  according  to  the  nature  of  the 
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Reynolds  number  effect,  scale  effects  should  not  be  expected  for  any 
of  the  first-order  derivatives,  e.g.,  Y^,  Y^,  Y^,  N^,  N^,  N^,  etc.,  which 
in  general  only  represent  lift  slope  characteristics.  In  the  ease  of  the 
higher  order  deriva wives,  however,  the  possibility  of  scale  effects  should 
be  considered,  as  it  is  likely  that  these  coefficients  would  be  influenced 
if  separation  or  flow  breakdown  occurred.  Normally,  higher  order  deriva¬ 
tives  of  the  motion  parameters  v  and  r,  for  instance  Y  ,  Y  etc.,  are 

vw  rrr 

determined  for  relatively  small  values  of  v  and  r  corresponding  to  angles 
of  attack  before  any  separation  effect  takes  place.  For  this  reason  scale 
effects  would  probably  be  negligible  also  for  these  coefficients.  This  is 
not  true  for  the  rudder,  as  the  rudder  deflection  for  which  rudder  character¬ 
istics  are  measured  also  will  cover  the  range  of  rudder  breakdown.  For  the 
derivatives  Y^gj  and  In  particular,  a  rational  correction  for  scale 

effects  should  be  considered. 

The  maximum  lift  is  sensitive  to  surface  roughness,  especially  near 
the  leading  edge.  Thus,  model  rudders  should  be  finished  as  smooth  as 
posjible  in  order  to  operate  in  a  well-defined  condition  and  to  obtain  re¬ 
peatable  measurements.  Similarly,  the  surface  roughness  of  the  full-scale 
rudder  should  be  taken  into  consideration  and  corrected  for  as  part  of 
the  above-mentioned  correction  of  rudder  derivatives  Y^^  and  for 

Reynolds  number  effect. 

Model  tests  should  be  carried  out  for  propeller  revolutions  corres¬ 
ponding  to  the  ship  propulsion  point  and  not  to  the  model  propulsion  point, 
which,  e.g.,  normally  would  have  to  be  applied  using  free-running,  self- 
propelled  models.  The  propeller  slipstream  can  thus  be  correctly  modeled. 
This  has  been  found  to  be  very  important,  as  it  has  a  great  effect  not 
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only  upon  the  rudder  derivatives  and  ,  but  also  upon  the  hull  deriva¬ 
tives  Yy,  Yr,  Ny  and  Nr* 

As  outlined  previously,  the  coefficients  X  ,  X  and  X  in  this 

VI  uu  uuu 

computer  program  are  calculated  on  the  basis  of  the  proper  ship  resistance 
values  and  a  power  assumption  corresponding  to  the  engine  setting  which 
would  be  attempted  during  an  actual  maneuver.  As  these  coefficients  are 
of  prime  importance  in  obtaining  the  correct  speed  reduction  during  a  maneu¬ 
ver,  it  is  found  that  a  principal  scale  effect  problem  has  thus  been  taken 
properly  into  account.  This  procedure  would  be  contrary  to  the  free-running 
model  techniques  where  the  difference  between  model  and  ship  resistance 
would  be  a  serious  problem  and  result  in  the  measurement  of  a  too  small 
speed  reduction  in  model  scale. 

The  foregoing  discussion  of  factors  influencing  scale  effect  should  in¬ 
dicate  that  It  is  possible  to  take  scale  effect  problems  into  account  in 
the  determination  of  the  different  coefficients  for  the  mathematical  model. 
Present  experience  might  be  insufficient  to  introduce  a  correction  for 
Reynolds  number  effect  as  suggested  for  the  rudder  derivatives  y666  and 
nevertheless,  a  correction  is  thought  to  be  feasible.  It  is  empha¬ 
sised  that  this  is  in  contrast  to  the  free-running  model  technique,  where 
the  scale  effect  problems  caused  by  incorrect  propulsion  point,  Reynolds 
number  effects,  etc.,  would  be  completely  mixed  up  in  the  model  results, 
leaving  only  very  little  room  for  introduction  of  scale  effect  corrections 
based  on  a  proper  physical  understanding  of  the  problem. 

VARIATIONS  OP  COEFFICIENTS  WITH  SPEED 

The  computer  program  has  been  based  on  a  solution  of  the  mathematical 
model  in  nondimensional  form;  consequently,  the  coefficients  used  as  input 
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data  to  the  program  should  be  applied  in  their  corresponding  nondimensional 
form. 

The  calculation  of  a  full-scale  trajectory  of  a  ship  maneuver  is  based 
on  d imansional izing  by  the  instantaneous  forward  velocity  u(t);  see  Equa¬ 
tions  (14)  and  (15).  When  a  certain  speed  loss  takes  place  during  a  maneuver, 
forces  and  moments  are  thus  basically  considered  as  being  proportional  with 

the  instantaneous  speed  squared,  and  coefficients  such  as  Y  ,  Y  ,  Y  , 
v  *  vu’  vuu’  ru’ 

Yruu>  Y^u>  Y^uu,  etc.,  which  represent  the  change  of  forces  and  moments  with 
speed,  should  only  reflect  the  extent  to  which  this  proportionality  does  not 
hold  true. 

Measurements  of  the  nondimensional  coefficients  Y  * ,  Y  1 ,  K  * ,  and  N  1 

v  r  v  r 

carried  out  for  various  ship  models  at  different  speed  values  have  indicated 


that  these  coefficients  are  largely  independent  of  speed.  Thus  coefficients 

vu  ’  vuu  ’  ru  *  ruu  ’  vu  *  vuu  *  ru  *  ruu  *  r 

the  change  with  speed,  are  negligible.  Consequently,  at  present  it  has  been 

found  reasonable  to  eliminate  these  coefficients  in  the  computer  program. 

For  the  rudder  derivatives  Y^’  and  1 ,  a  noteworthy  effect  has  been 
measured  for  a  change  in  forward  speed  especially  on  ships  where  the  rudder 
Is  situated  in  the  propeller  slipstream.  Apparently,  this  is  due  to  the 
fact  that  the  propeller  slipstream  is  nearly  constant  even  for  a  considerable 
change  of  forward  Speed,  because  propeller  revolutions  are  kept  more  or  less 
constant  during  a  maneuver.  Thus,  the  velocity  of  the  inflow  to  the  rudder 


Is  not  dependent  on  forward  speed  alone }  consequently,  the  nondimensional 
coefficients  Y^'  and  N^'  must  vary  as  a  function  of  forward  speed.  The 
coefficients  Ygu'  and  N«u' ,  which  represent  the  first  order  change  of  the 
rudder  derivatives  with  speed,  are  for  this  reason  thought  to  be  of 


considerable  importance,  and  they  should  be  included  in  an  experimental 
determination  of  the  various  coefficients. 


The  coefficients  Y^J,  representing  only  the  second-order  change 

of  Y^* and  N^1  with  speed,  have  nevertheless,  been  considered  negligible  and 
eliminated  in  the  program. 


The  coefficients  xrrB'»  X$5W’»  xv$u'>  a«d  X£ru'  i»  the  X-eqma- 

tion,  which  represent  the  change  of  X^',  X^',  Xfe',  X^/,  X^’,  and  X^'  with 
forward  speed,  have  similarly  been  omitted  from  the  computer  program  as 


they  are  thought  be  be  of  minor  importance  at  least  in  comparison  with  the 

dominating  coefficients  X  ,  X  ,  and  X 

u’  uu*  uuu 


RESUME  OF  COEFFICIENTS 


The  mathematical  model  developed  in  Equations  (10)  include  17  coeffi¬ 
cients  in  the  X-equation  and  24  coefficients  in  each  of  the  Y-  and  N-equa- 
tions.  As  mentioned  in  the  previous  section,  several  of  the  coefficients 
representing  change  of  nondimensional  forces  and  moments  with  forward  speed 
have  been  found  negligible  and  are  eliminated  in  the  computer  program. 

Obviously,  coefficients  are  of  varying  importance  with  respect  to  the 
accuracy  of  a  prediction,  and  a  classification  of  the  coefficients  has 
been  attempted  in  the  tunwary  of  the  coefficients  given  in  Tables  1-3, 
pages  25-27. 

The  tables  also  show  the  identifiers  that  have  been  used  for  the 
coefficients  in  the  computer  program  as  well  as  nondimensional  factors 
and  examples  of  the  numerical  values  taken  from  Reference  5.  The  planar 
motion  mechanism  test  technique,  which  could  be  used  to  measure  the  coeffi¬ 
cients,  is  mentioned  briefly. 
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Table  1  -  Summary  of  Coefficients  in  X-Equation 


■ . -"■-’I 

Variable 

X  -  fc 

q  >»  0  t  i  o  n 

— ■ - i 

Planar  Motion  Mechanism  Teat  Technique! 

Tovler 

Identifier 

For-ii’-. 

Relative 

or 

Calculation  method 

Expansion 

in 

lector 

Importance 

And  Cvr.amic 

FORI-Uh 

■h-UK 

of 

P~opor.se 

Frog ram 

Cooff. 

Tores 

Example 

(1) 

(E) 

(?) 

L _ 

>1 

<n-X.) 

X  UtOT 

I 

Estimated  from  theory  X.~-o.05  m 

Su 

X  ’j 

i  p  i.r-F?u 

-ko.q 

I 

1 

Calculated  on  the  baain  of  ship 

f  *uu 

x  ini 

i  p  Ltf? 

45.0 

I 

•  EKP-data  and  reaults  from  open- 

1  * 

7  *auu 

x  -a* 

i  f  i:rVu 

-10,5 

I 

water  propeller  to8t. 

v2 

Kv 

X  vv 

if  is/ 

-391.  3 

sex 

r 

tatlc  drift  ancle  teat 

r2 

(|  Xrr*m*c) 

X  RR 

iplEf4 

13.0 

«?i 

Pure  van  (angular  motion)  test  j 

4  2 

1*45 

X  ID 

2  7 

i  purV 

-94.3 

KI 

r 

tatic  drift  angle  teat 

O 

v  m 

1  x 

2  vvu 

j  ptEP2/u 

T  4U 

JL  t 

2  *rru 

i  p  LBP4/u 

0?AU 

1  y 

zh&u 

iptEP?u 

vr 

(  v* ) 

X  VR 

j  pLEF" 

7«3.0 

t! 

Taw  and  drift  angle  teat  -  m  is  known 

v6 

XV& 

X  VD 

i  pLEFCu 

N 

Static  drift  angle  teat 

r6 

Xr6 

X  RD 

i  plEP?u 

if.':'. 

Yaw  and  rudder  angle  test 

vrau 

Xvru 

i  pLBF-’/u 

HH 

viau 

Xv4u 

i  p  IBP2 

riau 

Xr6u 

i  pLEP3 

HB1 

- 

X* 

m 

7  7 
?  p  LBP  tr 

B 

Statir  drift  angle  teat 

- - - - 

(1)  The  Fortran  program  floes  not  include  all  terms  in  the  mathematical  model,  Equations  (10),  Certain 
coe  ficlents  have  been  left  out,  as  thov  have  been  considered  unimportant  for  the  accuracy  of  the 
pred ictlona. 

(?)  The  nondimenaional  coif ficients  have  been  taken  from  Reference  5. 

(1)  The  coefficients  have  been  divided  into  three  grades  according  to  their  importance  for  the 
accuracy  of  a  prediction.  The  most  important  coefficients  are  indicated  by  I;  coefficients  of 
minor  importance  by  Ml;  coefficients,  which  apparently  are  negligible,  by  N. 
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Table  2  -  Summary  of  Coefficients  in  Y-Equation 


Variable 

Y-Equatic 

n 

Planar  Motion  Mechanism  Test  Technique 

Taylor 

Identifier 

Nondlm. 

Nondlm. 

Relative 

Calculation  Method 

Expansion 

in 

Factor 

Coeff . 
»io5 

Importance 

And  Dynamic 

FORTRAN 

of 

Response 

Program 

from 

Cooff. 

Terms 

Example 

(1) 

(2) 

(3) 

MM 

(m-Y.) 

Y  VDOT 

ip  LBP3 

I 

Pure  sway  (transverse  motion)  test 

n 

<»vV 

Y  ROOT 

J  p  LBP4 

§|§j 

I 

Pure  yaw  (angular  motion)  test 

V 

*v 

Y  V 

i p  LBF1 2 3u 

-1160.4 

I 

Static  drift  angle  test 

v3 

Kvv 

Y  VW 

i  p  LBP2/u 

-S07B.2 

MI 

Static  drift  angle  teat 

vr2 

1  Y 

2  vrr 

Y  VRR 

i  p  LBP4/u 

0.0 

N 

Yaw  and  drift  angle  test 

vA2 

1  Y 

2  v66 

Y  VDD 

i p  LBP2u 

-3.B 

N 

Static  drift  angle  test 

VAU 

i  p  IBP2 

2 

tau 

—  Y 

2  «uu 

_ 

J  p LBP2/u 

r 

(  Y  -mu) 
r 

i f LBP3u 

-499.0 

i 

Pure  yaw  (angular  motion)  test 

r3 

1  v 

7>  rrr 

i p IBP5/u 

0.0 

N 

Pure  yaw  (angular  motion)  test 

rv2 

•1  Y 

2  rw 

i p LBP3/u 

15356.0 

I 

Yaw  and  drift  angle  test 

r62 

KaA 

B 

i p  LBP3u 

0,0 

N 

Yaw  and  rudder  angle  test 

rau 

Y 

ru 

ip  IBP3 

2 

rau 

A  y 

2  ruu 

■ 

i  p  IBP3/u 

A 

*6 

Y  D 

i  p  LBP2u2 

277.9 

I 

Static  drift  angle  test 

6  3 

\  *466 

Y  DDD 

i  p  LBP2u2 

■EH 

MI 

Static  drift  angle  test 

6v2 

1  Y 

2  Aw 

Y  DVV 

i  p  EBP2 

MI 

Static  drift  angle  test 

Ar2 

1  *6rr 

Y  DRR 

i  p  LBP4 

Hll 

a 

Yaw  and  madder  angle  test 

Aau 

*6u 

Y  DO 

i  p LBP2u 

(OeO) 

MI 

Static  drift  angle  teet 

executed  at  varloue  speed  volues 

Aau2 

1  y 

2  Auu 

i  p  LBP2 

vrfc 

*vrA 

Y  VRD 

i  pLBP3 

N 

Yaw  and  drift  angle  teat 

executed  at  various  speed  values 

• 

** 

I. 

i  f  LBP2u2 

-3.6 

MI 

Static  drift  angle  test 

AU 

Y*u 

1 

i  p  lbp2u 

(0.0) 

N 

Static  drift  angle  test 

executed  at  various  speed  values 

AU2 

V*uu 

EM 

i  p  LBP2 

(1)  The  FORTRAN  proprom  does  not  incorporate  all  terms  in  the  mathematical  model,  Equations  (10), 
Certain  coefficients  have  been  left  out,  as  they  have  been  considered  without  importance  for  the 
accuracy  of  the  predictions. 

(2)  The  nondimenaional  coefficients  have  been  taken  from  Reference  5  except  values  enclosed  in 
parenthesis,  for  which  nn  data  were  available. 

(3)  The  coefficients  have  been  divided  into  three  grades  according  to  their  importance  for  the 
accuracy  of  a  prediction.  The  moat  important  coefficients,  which  should  be  available  in  order 
to  obtain  a  prediction,  are  marked  by  Is  coefficients  of  minor  importance  by  Mis  coefficients 
which  apparently  are  negligible,  by  N. 
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Table  3  -  Summary  of  Coefficients  in  N-Equation 


Variable 

N  -  E 

q  u  a  t  i 

n 

Planar  Motion  Mechanism  Test  Technique 

Taylor 
Expansion 
And  Dynamic 
Response 
Terms 

Identifier 

in 

FORTRAN 

Program 

Nondim. 

Factor 

Nondiv* 
Coeff • 
«105 
from 
Example 

Relative 

Importance 

of 

Coeff. 

or 

Calculation  Method 

(1) 

(?) 

(3) 

V 

("vV 

N  VDOT 

J  p  LBP4 

-22.7 

1 

Pure  swav  (transverse  motion)  test 

r 

(I  -N. ^ 

Z  I 

K  RDOT 

4  9  LB?5 

82.9 

I 

Pure  >-aw  (angular  motion)  test 

V 

N 

V 

N  V 

4  p  lbp1 * 3u 

-263.5 

I 

Static  drift  angle  teat 

v5 

I  s 

6  vvv 

N  VVV 

4  p  LBP5/ u 

1636.1 

MI 

Static  drift  angle  test 

vr£ 

Ik 

2  vrr 

K  VRR 

4  9  lbp5/u 

0.0 

N 

Yaw  and  drift  angle  test 

v62 

1  ,, 

2 

N  VDD 

4  9  LBP3u 

12.5 

N 

Static  drift  angle  test 

V.1U 

Kvu 

4  p  LBP1 

2 

vi»u 

I  Hvuu 

4  9  lbp3/u 

r 

(  Kr-mx0u) 

N  R 

4  9  LBP4u 

-166.0 

i 

Pure  yaw  (angular  motion)  test 

r1 

l  Krrr 

K  RKR 

4  p  LBP6/u 

0.0 

» 

Pure  vaw  (angular  motion)  te3t 

rv? 

1  Hrvv 

N  RW 

4  p  LBP4/ u 

-5483.0 

; 

Yaw  and  drift  angle  test 

r{2 

?"r44 

H  RDD 

4  p  lbp4u 

0  0 

N 

Yaw  and  rudder  an*le  test 

r«u 

"ru 

4  9  LBP4 

rau2 

1  N 

2  ruu 

4  9  lbp4/u 

6 

k6 

H  D 

4  9  lbpV 

-138.8 

I 

Static  drift  angle  test 

41 

Z  86to 

N  DDD 

4  p  lbp  V 

45.0 

MI 

Static  drift  angle  test 

6v2 

f  M6w 

H  DVV 

4  9  LBP3 

-489.0 

MI 

Static  drift  angle  tost 

6r? 

lN6rr 

H  DRR 

4  P  LEP4u 

n.O 

B 

Yaw  and  rudder  angle  test 

6au 

H6u 

K  DU 

4  9  lep3u 

(0.0) 

MI 

Static  drift  angle  test 

executed  at  various  speed  values 

6au2 

2  !'4uu 

4  9  LBP3 

vri 

Kvri 

H  VRD 

4  9  LBP4 

0.0 

B 

Yaw  and  drift  angle  test 

executed  for  various  speed  values 

- 

B  0 

4  9  LBP3u2 

2.9 

MI 

Static  drift  angle  test 

AU 

N  OU 

4  9  LEP?u 

(0.0) 

B 

Static  drift  angle  test 

executed  at  various  speed  values 

AU2 

*uu 

4  p  LBP3 

(1)  The  FORTRAN  program  doen  not  incorporate  oil  tetrns  in  the  mathematical  model,  Equations  (10). 
Certoin  coefficients  have  teen  left  out,  as  they  have  been  considered  without  importance  for  the 
accuracy  of  the  predictions. 

(?)  The  nondimensional  coefficients  have  been  taken  from  Refererce  5  except  values  enclosed  in 
parenthesis,  for  which  no  data  were  available. 

(3)  Tho  coefficientn  have  been  divided  into  three  rrndec  according  to  their  importance  for  the 
accuracy  of  a  prediction.  The  most  important  coefficients,  which  should  be  available  in  order 
to  obtain  a  prediction,  are  marked  by  X;  coefficients  of  minor  importance  by  Kl;  coefficients, 
which  apparently  are  negligible,  tv  N. 
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COMPUTER  PROGRAM  FOR  PREDICTION  OF  STANDARD  MANEUVERS 


The  solution  of  the  mathematical  model  for  steering  and  maneuvering 
has  been  programmed  in  the  FORTRAN  II  language  available  for  the  IBM  7090 
computer  at  '3MB.  The  program  is  designated  AML  Problem  XPMC.  The  FORTRAN 
listing  of  the  computer  program  is  included  in  Appendix  C  of  this  report. 

INPUT  DATA 

Data  forms  have  been  worked  out  to  help  in  the  accurate  preparation  of 
input  data  for  the  computer  program.  An  example  of  the  data  forms  is  given 
in  Appendix  A,  and  the  following  discussion  of  the  input  data  refers  to  Chi' 
example . 

The  input  data  consist  of  two  parts;  (1)  Specification  data,  page  52 
find  (2)  Ship  data,  pages  53-56  . 

Specification  Data 

The  specification  data  describe  the  maneuvers  which  should  be  predictec 
at  the  execution  of  the  program.  Four  different  types  of  calculations  can 
be  specified  and  carried  out  by  the  program; 

1.  Calculation  of  the  turning  circle  parameters  as  defined  in  Figure 
2.  The  parameters  are  calculated  for  a  series  of  different  rudder 
deflections,  which  should  be  specified  in  the  data  form. 

2.  Calculation  of  the  turning  circle  trajectory  for  a  certain  rudder 
deflection.  Parameters  such  as  advance,  transfer,  speed,  heading 
angle,  angular  velocity,  and  drift  angle  are  presented  on  a  time 
basis  for  each  10  sec  until  a  A50-deg  turn  has  been  executed.  The 
turning  circle  calculation  can  be  specified  for  several  rudder 
deflections  at  each  execution  of  the  program. 
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3.  Calculation  of  the  zig-zag  maneuver  as  defined  diagrammc',:ically 
in  Figure  3.  The  same  parameters  as  mentioned  above  for  the 
turning  circle  calculation  are  presented  on  the  basis  of  a  time 
interval  of  10  sec.  The  calculation  of  the  zig-zag  maneuver  can 
be  repeated  for  different  limits  of  the  rudder  and  heading  angle 
at  each  execution  of  the  program  if  this  is  desired. 

4.  Calculation  of  spiral  maneuver.  This  maneuver  is  executed  as 
usual  starting  with  a  specified  positive  rudder  deflection,  step¬ 
wise  reducing  the  rudder  angle  to  a  specified  negative  rudder 
deflection  and  vice  versa.  To  obtain  an  accurate  determination 
of  a  possible  Loop  phenomenon,  a  smaller  difference  between  con¬ 
secutive  rudder  positions  can  be  specified  in  the  range  around 
zero  rudder  deflection  (see  figure  on  data  form,  page52). 

Port  or  left  rudder  is  considered  a  positive  rudder  deflection  in 
the  program.  Similarly,  starboard  or  right  rudder  corresponds  to  a  nega¬ 
tive  deflection.  Hie  rudder  deflections  should  be  specified  accordingly 
in  the  data  forms. 

The  35-deg  turning  circle,  20-20  deg  zig-zag,  and  spiral  maneuvers 
are  referred  to  as  the  "Standard  Maneuvers"  which  are  used  to  evaluate 
performance  qualities  of  a  surface  ship.  The  maneuvers,  which  have  been 
specified  on  the  example  of  the  data  form  in  Appendix  A,  actually  corres¬ 
pond  to  those  "Standard  Maneuvers." 

A  graphic  display  of  the  computer  results,  that  is,  turning  circle 
trajectory,  zig-zag,  and  spiral,  maneuver  can  be  obtained  directly  from 
the  TMB  computer  by  means  of  the  on-line  Charactron  plotting  equipment. 
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Such  a  plotting  of  the  results  can  be  specified  on  the  data  form  as  a  pert 
of  the  specification  data. 

Ship  Data 

The  ship  data  have  been  divided  into  three  groups:  (1)  Principle 
ship  data,  page  53,  (2)  EHP-data  and  open- water  propeller  characteristics, 
pages  54-55  ,  and  (3)  nondimens ional  coefficients,  page  56. 

The  principle  ship  data  include  particulars  such  as  ship  length,  beam, 
draft,  displacement,  propeller  dimensions,  wake  coefficient,  thrust  deduc¬ 
tion  coefficient,  etc.  This  group  of  data,  furthermore,  incorporates 
values  for  the  rudder  system,  such  as  rudder  rate  and  timelag  discussed 
previously  in  the  section  "Definition  of  Rudder  Deflection,"  page  16. 

The  ship  effective  horsepower  data  and  open-water  characteristics 
for  the  propeller  (Data  Group  2)  together  with  data  for  approach  speed, 
wake  coefficient,  thrust  deduction  coefficient,  and  information  about  the 
type  of  the  ship  propulsion  plant '(Data  Group  1)  are  the  basis  for  calcu¬ 
lating  the  coefficients  X  ,  X  .  and  X  .  ,  as  discussed  in  the  section 
"Calculation  of  Coefficients  in  X-Equation,"  page  19.  The  ship  effective 
horsepower  data  should  be  given  for  a  range  of  speed  values  covering  the 
values  to  be  encountered  during  the  maneuvers.  The  roughness  or  extra¬ 
polation  allowance  used  in  the  preparation  of  these  data  should  correspond 
to  the  condition  of  the  ship  hull  roughness  for  which  the  maneuvering  pre¬ 
dictions  are  desired.  Similarly,  the  wake  and  thrust  deduction  coeffi¬ 
cients  should  be  actual  ship  values;  for  Instance,  those  obtained  from 
model  tests  and  corrected  for  possible  scale  effects. 

The  calculation  of  the  Xu,  Xuu,  and  Xyuu  coefficients  can  be  based 
either  on  the  assumption  that  propeller  revolutions  will  be  kept  constant 


30 


during  the  maneuvers  or  on  the  assumption  that  propeller  torque  will  vary 
proportionally  to  propeller  revolutions  in  a  certain  power.  In  the  first 
mentioned  case,  the  input  value  named  TYPE  on  the  data  form  should  be 
chosen  as  a  value  smaller  than  -5.0,  whereas  in  other  cases  TYPE  represent 
the  power  factor  in  the  propeller  torque  equation: 

V\  -  <nt/ntl)TyPE  <23> 

The  proper  value  of  the  power  factor  depends  upon  the  type  of  power  plant 
and  should  be  estimated  from  information  about  the  actual  ship  under  con¬ 
sideration.  If  TYPE=0,  Equation  (23)  would  represent  the  case  where  torque 

remains  constant  during  a  maneuver  since  Q  =  Q  .  This  would  largely 

C1 

correspond  to  the  conditions  of  a  Diesel  engine,  which  would  develop  a 
constant  torque  independent  of  speed  and  propeller  revolutions  so  long 
as  the  engine  setting  is  kept  constant.  If  TYFE=-1.0,  Equation  (23)  would 
change  to  Qt*nt=  Qt  *nt  »  which  actually  would  correspond  to  a  turbine  cap¬ 
able  of  maintaining  a  constant  power  output  (which,  for  instance,  would  be 
the  case  for  a  turbine  with  semi-automatic  throttle  control) . 

The  nondimensional  coefficients  to  be  stated  as  the  last  group  of  the 
ship  data  (Data  Group  3)  are  those  coefficients  obtained  from  captive 
model  testing,  previously  discussed  in  the  Chapter  "Coefficients  in  Mathe¬ 
matical  Model,"  page  17.  Tables  1,  2,  and  3,  pages  25-27,  show  the  re- 
lationship  between  the  hydrodynamic  derivatives  in  the  mathematical  model, 
Equation  (10) ,  and  the  corresponding  identifiers  used  in  the  program,  and 
on  the  data  form.  The  tables  also  give  the  nondimensionalizing  factors, 

S* 

to  be  used  in  the  data  preparation. 


31 


the  data  forms  in  Appendix  A  give,  as  an  example,  values  of  the  various 


ship  data  corresponding  to  a  cargo  chip*  Except  for  coefficients  Y^u  and 
N^u,  these  data  have  been  taken  from  the  data  and  coefficients  published  in 
Reference  5  for  the  MARINER  hull  form.  The  coefficients  Y^u  and  N^u  are 
thought  to  be  of  some  importance  for  the  accuracy  of  the  predictions,  but 
they  have  been  put  equal  to  zero  in  the  data  forms  as  no  model  tests  are 
available  at  present. 

The  data  forms  are,  in  general,  thought  to  be  self-explanatory.  They 
contain  the  FORMAT  specifications  which  necessarily  must  be  known  for  the 
preparation  of  the  punched  cards. 


OUTPUT  FORM 

The  output  from  the  computer  program  is  presented  in  the  form  of  a 
printed  "prediction  report"  and,  if  desired,  as  graphs  plotted  by  means  of 
the  Charactron  Microfilm  Recorder.  An  example  of  both  types  of  output  is 
given  in  Appendix  B,  which  presents  the  results  corresponding  to  the  input 
data  shown  in  the  data  forms.  Appendix  A.  The  following  discussion  of  the 
output  refers  to  the  example  in  Appendix  B. 

Prediction  Report 

The  first  pages  of  the  prediction  report  define  precisely  the  input 
data  on  the  basis  of  which  the  prediction  has  been  carried  out.  On  PAGE  1 

it  gives  the  principal  ship  data;  on  PAGE  2,  the  EHP-data  and  open-water 

% 

propeller  curves;  and  on  PAGE  3,  the  nondimensional  hydrodynamic  coeffi¬ 
cients;  see  pages  58-60. 

PAGE  2  of  the  output  shows  the  calculation  of  the  coefficients  Xu> 

Xuu,  and  Xuuu,  which  in  this  case  has  been  carried  out  under  the  assumption 


of  a  constant  power  output  from  the  turbine  (this  means  that  propeller 
torque  multiplied  by  propeller  revolutions  has  been  kept  constant  for  the 
different  values  of  speed) .  It  is  seen  that  this  power  assumption  gives  a 
slight  variation  of  the  propeller  revolutions,  varying  from  68.6  rpm  at  the 
15-knot  approach  speed  to  55.5  at  a  speed  of  7.0  knots. 

To  facilitate  a  straightforward  evaluation  of  the  inherent  dynamic 
stability  of  the  ship,  the  output  includes  on  PAGE  3  values  for  the  non- 
dimensional  stability  criterion  and  stability  roots  as  well  as  the  slope 
of  the  r-6  curve  in  sec"’*'.  These  quantities  have  been  computed  on  the  basis 
of  the  linear  theory  according  to  which  the  criteria  for  dynamic  stability 
as  mentioned  in  Equation  (9)  become: 


C  =  Yv(N.-mxGu1)  -  Nv(Yr-mUl)  >  0  (9) 

The  stability  roots,  which  all  should  be  negative  for  a  stable  ship,  are, 
in  accordance  with  Reference  3,  defined  by 


ITji  f  SIGMA  1^ 
Suj  =  I  SIGMA  2  J  “ 
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where : 


A  =  (m-Y^)  (I#-N.)  -  (taxG»Yf)  (mXg-N.)  , 

B  =  -(m-Y^)(»r-m*Gu1)  -  <IJ,-N^)Yv+  (mxG~Y^)Nv+  faXg-JL  )  (Y^mUj) ,  and 
C  =  stability  criterion.  Equation  (9) . 

The  slope  of  the  r-6  curve,  which  represents  the  change  in  angular 
velocity  r  for  a  small  rudder  deflection  S,  is  expressed  by 
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The  slope  is  negative  for  a  stable  ship,  infinite  for  the  marginally  stable 
ship,  and  positive  for  the  unstable  ship  as  indicated  in  Figure  4. 

The  next  pages  of  the  prediction  report  give  the  results  fron  the  four 
different  types  of  calculation  which,  as  described  previously,  can  be 
carried  out  by  the  program: 

1.  Calculation  of  turning  circle  parameters,  PAGE  4 

2.  Calculation  of  turning  circle  for  specified  rudder  angle,  PAGE  5 

3.  Calculation  of  zig-zag  maneuver,  PAGE  7 

4.  Calculation  of  spiral  maneuver,  PAGE  9. 

It  should  be  noted  that  the  speed  values  given  in  the  results  correspond 

-*>  H? - T 

to  the  velocity  vector  U  =  uu“  +  v  and  not  to  the  forward  component 
u  of  the  velocity. 

Charactron  Microfilm  Plotting  of  Maneuvers 

The  Ch&ractron  Plotting  Equipment,  which  is  available  as  an  on-line 
output  facility  at  1MB,  permits  the  plotting  and  recording  of  results 
directly  on  microfilm,  to  be  developed  and  enlarged  subsequently.  Pages 
68-70,  Appendix  B,  show  examples  of  the  three  types  of  plots  that  can  be 
obtained  from  the  computer  program  in  this  way. 

The  plot  of  the  turning  circle  trajectory,  page  68,  gives  the  path 
of  the  origin  of  the  ship.  The  origin  of  the  coordinate  axis  system 
corresponds  to  the  point  of  rudder  execute.  Hie  plotting  is  continued 
Until  a  540-deg  turning  circle  has  been  completed. 

The  plot  of  the  zig-zag  maneuver,  page  69,  gives  the  well-known  pre¬ 
sentation  of  rudder  angle  and  heading  angle  on  a  time  basis.  Rudder 
positions  are  indicated  by  an  asterisk  plotted  every  10  sec.  Similarly, 
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the  heading  angle  is  plotted  with  10-sec  time  intervals,  but  in  this  case 
a  straight  line  connecting  subsequent  points  produces  the  "continuous" 
curve  seen  in  the  graph.  In  this  connection,  it  should  be  noted  that  the 
program  is  based  on  a  solution  of  the  mathematical  model  using  a  time 
interval  of  1  sec,  as  mentioned  on  page  16.  However,  the  plotting  of  the 
maneuvers  has  been  based  on  points  with  time  intervals  of  10  sec  only  in 
order  to  reduce  the  amount  of  data  to  be  stored  in  the  computer.  The 
somewhat  stepwise  appearance  of  the  heading  angle  curve  is,  for  this  reason, 
due  to  the  method  of  plotting,  and  should  not  be  taken  as  an  expression  for 
the  accuracy  of  the  computer  solution. 

The  result  from  the  spiral  maneuver,  page  70 ,  is  presented  as  a  plot 
of  rate  of  change  of  heading  in  degrees  per  second  versus  the  different 
rudder  positions.  The  results  are  plotted  as  discrete  points  only,  and  it 
might  be  necessary  to  consult  the  printed  results  in  order  to  separate 
points  obtained  during  the  spiral  maneuver  for  decreasing  and  increasing 
rudder  angle,  respectively.  (The  line  connecting  the  points  on  the  example 
has  not  been  drawn  by  the  recorder,  but  has  been  inserted  afterward  by  hand 
to  help  in  reading  the  points  on  the  graph.) 

RESULTS  OF  SAMPLE  CALCULATIONS 

The  output  example  presented  In  Appendix  B  has  been  calculated  on  the 
basis  of  hydrodynamic  coefficients  for  the  MARINER  form  measured  by  planar 
motion  mechanism  tests  and  reported  in  Reference  5.  The  comparison  between 
full-scale  trials  and  computer  predictions  given  in  Reference  5  shows  that 
the  computer  solution  of  the  nonlinear  mathematical  model  presents  an  ac¬ 
curate  method  for  the  prediction  of  the  "Standard  Maneuvers."  No  measure- 


Bents  were  available  for  the  coefficients  Y^u  and  N^u,  which  represent  the 
first-order  change  of  rudder  derivatives  and  with  speed.  It  is 

however,  that  inclusion  of  values  for  these  coefficients  would 
have  improved  the  accuracy  further,  particularly  in  the  prediction  of 
tight  maneuvers  where  a  considerable  speed  loss  takes  place. 

To  demonstrate  the  potential  of  the  computer  program  and  its  ability 
to  give  detailed  information  about  the  maneuvers,  different  sample  calcu¬ 
lations  have  been  worked  out  and  are  presented  in  the  following  sections. 

INFLUENCE  OF  TYPE  OF  POWER  PLANT  UPON  THE  SPEED  LOSS  IN  MANEUVERS 

The  type  of  power  plant  has  a  considerable  effect  upon  the  speed  loss 
which  takes  place  during  a  maneuver.  To  show  this  influence,  predictions 
of  the  "Standard  Maneuvers"  have  been  calculated  for  three  different  cases 
assuming. constant  propeller  revolutions,  constant  engine  power  (turbine 
ship),  and  constant  torque  (Diesel  ship),  respectively.  All  three  set6 
of  predictions  have  been  carried  out  on  the  basis  of  the  MARINER  coeffi¬ 
cients  for  an  approach  speed  of  15  knots.  Thus,  the  prediction  for  the 
turbine  ship  corresponds  to  the  results  of  the  example  given  in  Appendix  B. 

Figures  5  through  8  present  some  of  the  results  obtained  from  the 
three  predictions.  Figure  5  shows  change  of  propeller  revolutions  as  a 
function  of  forward  velocity.  Figure  6  gives  velocity  turn  entry  trans¬ 
ient  for  the  35-deg  port  rudder  turning  circle,  and  Figure  7  gives,  simi¬ 
larly,  change  in  velocity  predicted  for  the  zig-zag  maneuver.  In  Figure 
8,  the  results  from  the  spiral  maneuver  have  been  presented  in  the  usual 
form  as  rate  of  change  of  heading  versus  rudder  angle. 
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The  figures  clear ly  indicate  that  the  speed  loss  is  greatly  influenced 
by  the  power  assumption.  For  the  Diesel  ship,  it  takes  considerably  more 
time  before  steady  conditions  are  obtained,  because  the  speed  loss  is  also 
considerably  greater  (Figures  6  and  7) . 

Free-running  model  tests,  which  often  would  be  carried  out  for  con¬ 
stant  propeller  revolutions,  would  apparently  indicate  a  smaller  speed  loss 
ftd  a  shorter  transition  period  than  full-scale  trials. 

The  trajectories  of  the  predicted  maneuvers  have  been  found  to  be 

independent  of  the  power  assumption.  This  would  not  have  been  the  case 

had  coefficients  for  Yt  and  Nt  been  included  in  the  set  of  coefficients 

ou  Ou 

used  for  the  predictions.  Nevertheless,  it  indicates  that  trajectories, 
in  general,  would  be  independent  of  the  speed  loss  encountered  during  a 
maneuver.  This  further  indicates  that  it  might  be  advantageous  to  compare 
and  evaluate  maneuvering  performance  on  the  basis  of  measurements,  which 
are  independent  of  the  speed  loss  and  consequently  are  independent  of  .-*• 
the  power  plant  in  the  ship.  Results  from  the  spiral  could  be  presented 
in  a  form  suggested  in  Figure  9  as  a  plot  of  the  reciprocal  of  the  turning 
radius  versus  rudder  deflection  instead  of  in  the  usual  graph  shown  in 
Figure  8,  which  is  influenced  by  the  power  assumption.  This  representation 
would,  in  general,  be  Independent  of  the  power  assumption,  which  might  be 
difficult  to  obtain  correctly  from  full-scale  trials.  Results  from  the 
three  sets  of  predictions  would  in  this  way  he  plotted  as  a  single  curve 
independent  of  the  speed  loss. 

Similarly,  evaluation  of  the  xig-zag  maneuver  on  the  basis  of  ''period” 
and  "reach"  (gee  Figure  3)  would  be  influenced  (but  only  slightly)  by  the 
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power  assumption.  An  evaluation  based  on  factors  independent  of  time  would 
be  preferable. 

PREDICTION  OP  ZIG-ZAG  MANEUVERS  FOR  DIFFERENT  VALUES  OF  SHIP  INERTIA. 

The  predictions  present ,d  in  Appendix  B  have  been  computed  on  the 

basis  of  an  approximate  value  for  the  ship  moment  of  inertia  I  .  The  non- 

z 

dimensional  value  for. the  inertia  has  been  taken  as  I  ' »  39.2  •  IQ-'’,  and 
the  nondimensional  coefficient  as  (IJ-NJ)  =  N  RDOT  =  82,9  -  10 .  Any  in¬ 
accuracy  in  this  value  would  in  particular  have  an  effect  on  the  prediction 
of  the  zig-zag  maneuver.  To  estimate  this  influence,  supplementary  calcu¬ 
lations  have  been  carried  out,  assuming  the  inertia  to  be  25  percent  larger 
and  smaller,  respectively.  The  effect  of  this  change  is  shown  in  Figure 
10.  The  characteristic  measures,  overshoot,  reach,  and  period  as  defined 
in  Figure  3,  are  influenced,  but  nevertheless  it  is  comforting  to  see  that 
even  a  considerable  error  in  the  estimation  of  ship  inertia  would  introduce 
only  a  small  change  in  the  maneuvering  qualities  of  the  ship. 

This  example  at  the  same  time  indicates  the  flexibility  of  the  predic¬ 
tion  method.  Model  testing  can  be  executed  for  any  value  of  model  inertia, 
because  the  appropriate  ship  value  can  be  introduced  at  the  time  of  data 
preparation  for  the  computer  program.  This  is  in  contrast  to  the  free- 
running  model  technique,  where  model  inertia  should  be  properly  scaled. 
Furthermore,  strictly  speaking,  model  results  would  correspond  to  only 
one  value  of  ship  inertia. 

PREDICTION  OF  LOOP  PHENOMENON  IN  SPIRAL  MANEUVER 

The  application  of  the  nonlinear  mathematical  model  makes  it  possible 
to  give  realistic  maneuvering  predictions  even  for  ships  which  are  dynami- 
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sally  unstable  on  a  straight  course.  This  is  illustrated  by  Figure  11, 
which  shows  results  from  the  spiral  maneuver  predicted  for  four  different 
hull  forms,  two  of  which  have  been  unstable  v'.ile  the  others  have  been 
marginally  stable  and  stable,  respectively.  The  figure  shows  plots  of 
the  rate  oi  change  of  heading  versus  rudder  angle,  and  it  is  seen  that 
the  unstable  ships  exibit  a  zone  in  which  there  is  a  lack  of  preferential 
rate  of  change  of  heading  with  rudder  angle.  The  "loop"  phenomenon  asso¬ 
ciated  with  an  unstable  hull  form  has  thus  been  reconstructed  exactly  by 
the  computer  program. 


Table  4  -  Nondimensional  Coefficients  Governing  the  Criteria  for 
Dynamic  Stability  for  Stable,  Marginally  Stable,  and 
Unstable  Hull  Forms 


Nondim. 

Coeff . 

Hull  Form 

Stable 

MARINER 

Hull  Form 

Marginally 

Stable 

Hull  Form 
Unstable 

Hull  Form 
Unstable 

Y  -105 

V  f- 

-1160.4 

-1044.0 

-928.0 

-812.0 

V10, 

-263.5 

-290.0 

-316.0 

-343.0 

V“ 

298.0 

268,0 

238.0 

209.0 

Nr-105 

-184.3 

-166.3 

-147.3 

-129.3 

O^-mu^.fc)5 

-499.0 

-529.0 

-559,0 

-588.0 

CNr-mxGUi)  •  105 

-166.0 

-148,0 

-129.0 

-111.0 

c-io5 

0.61 

0.01 

-0.57 

-1.12 

The  curve  representing  the  stable  ship  corresponds  to  the  results  pre¬ 
sented  in  previous  examples  predicted  on  the  basis  of  the  hydrodynamic  coef¬ 
ficients  for  the  MARINER  form.  The  results  for  the  marginally  stable  and 
unstable  hull  forms  have  been  obtained  on  the  basis  of  hydrodynamic  coeffi¬ 
cients  derived  from  the  MARINER  values  by  changing  the  four  coefficients 
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T  s  H  ,  ?  ,  Arid  N  ,  which  govern  the  criterion  for  dynamic  ethB^ity 
given  in  Equation  (9).  Thus,  in  order  that  C  becomes  zero  for  the  mar¬ 
ginally  stable  and  negative  for  the  unstable,  the  derivatives  have  been 
changed  10,  20,  and  30  percent  as  shown  in  Table  4* 

SLOPED  LOOP  PHENOMENON  IN  SPIRAL  MANEUVER 

It  is  important  in  the  execution  of  a  full-scale  spiral  maneuver  to 
wait  a  sufficient  period  of  time  until  steady  conditions  have  been  reached 
before  measuring  rate- of  change jof  heading,  speed,  etc.,  and  before  order¬ 
ing  the  next  rudder  deflection.  For  certain  rudder  positions,  it  might, 
however,  take  a  considerable  time  before  the  motion  becomes  steady  or  it 
might  be  difficult  to  recognize  that  the  ship  actually  is  in  a  transition 
period.  For  this  reason,  measurements  might  be  taken  too  hastily.  As  a 
result,  the  spiral  maneuver  can  exibit  a  sloped  loop  phenomenon  even  for 
a  ship  that  is  in  reality  perfectly  stable. 

To  illustrate  the  sloped  loop  phenomenon,  predictions  have  been 
carried  out  for  the  (stable)  MARINER  form  used  in  previous  examples, 
executing  the  spiral  maneuver  with  a  limited  time  interval  between  consec¬ 
utive  rudder  deflections.  The  spiral  maneuver  has  been  computed  in  two 
cases  using  time  intervals  of  60  and  120  sec,  respectively.  The  results 
from  these  predictions  are  shown  in  Figure  12  together  with  the  results 
from  the  spiral  maneuver,  where  no  time  limit  has  been  applied.  It  is 
seen  that  a  double  curve  or  sloped  loop  is  obtained  in  the  cases  where 
premature  measurements  have  been  taken. 

She  possibility  of  a  sloped  loop  phenomenon  should  be  kept  in  mind, 
especially  when  evaluating  the  full-scale  spiral  maneuver  results  from  a 
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ship  which  might  be  marginally  stable.  A  sloped  loop  would,  in  such  a 
case,  easily  be  interpreted  as  the  loop  associated  with  a  dynamically 
unstable  ship. 


CONCLUSIONS 

The  computer  program  permits  the  calculation  of  steering  and  maneu¬ 
vering  trials  of  surface  ships  giving  predictions  of  the  turning  circle, 
zig-zag  and  spiral  maneuvers.  Predictions  are  presented  in  the  form  of  a 
printed  report  and  graphs  plotted  by  the  on-line  Charactron  Microfilm 
Recorder , 

The  program  is  based  on  the  solution  of  a  nonlinear  mathematical 
model  describing  the  motion  of  a  ship  in  the  horizontal  plane.  The  mathe¬ 
matical  model  is  developed  from  the  equationr  of  motion  using  a  third- 
order  Taylor  expansion  of  forces  and  moments.  The  model  has  been  reduced 
to  a  solvable  form  on  the  basis  of  the  following  assumptions: 

1.  Influence  from  rolling  of  the  ship  is  negligible  upon  maneuvering 
predictions  (page  6)  . 

2.  Forces  and  moments  can  be  considered  to  be  symmetrical  except  for 
side  force  from  propeller  (pages  10  and  11) . 

3.  No  second  or  higher  order  acceleration  terms  can  be  expected  in 
the  Taylor  expansion  of  forces  and  moments.  Similarly,  cross¬ 
coupling  between  acceleration  and  velocity  parameters  is  negli¬ 
gible  (page  11) . 

4.  Change  of  nondimen sional  coefficients  Y^',  Y^,  Nv',  with  speed 
is  negligible  (page  23 )  . 

Input  to  the  program  can  be  prepared  by  means  of  data  forms.  The 
data  consist  of  the  hydrodynamic  force  and  moment  coefficients  measured  by 
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captive  model  test  technique,  ship  EHP-data,  open-water  propeller 
characteristics,  as  well  as  data  for  the  rudder  system  and  types  of  power 
plant  of  the  ship. 

The  combination  of  captive  model  testing  and  the  computer  prediction 
of  maneuvers  permits  scaling  laws  to  be  taken  into  account  in  a  proper 
fashion.  The  hydrodynamic  coefficients  can  be  obtained  from  captive  model 
tests  executed  with  the  model  propelled  at  the  ship  propulsion  point.  It 
is  further  emphasized  that  coefficients,  if  experience  is  available,  can 
be  corrected  for  Reynolds  number  effects.  Coefficients,  which  are  of 
principal  importance  for  the  determination  of  speed  loss  during  maneuvers, 
are  computed  in  the  program  from  ship  resistance  values^  eliminating  the 
scale  effect  problem  arising  because,  of  the  difference  between  ship  and 
model  resistance. 

The  application  of  the  nonlinear  mathematical  model  makes  it  feasible 
to  give  accurate  predictions  for  any  type  of  maneuvers  including  tight 
maneuvers.  Also  the  loop  phenomenon  associated  with  the  spiral  maneuver 
for  a  ship  which  is  dynamically  unstable  on  a  straight  course  is  readily 
obtained  by  the  program 


RECOMMENDATIONS 

Most  of  the  assumptions  on  which  the  computer  program  is  based  have 
been  shown  to  hold  true,  e.g.,  in  the  experimental  measurements  reported 
in  Reference  5.  However,  further  model  tests  are  recommended  in  order 
to  confirm  the  assumptions. 

Comparison  between  full-scale  trials  and  computer  predictions  given 
in  Reference  5  shows  a  promising  agreement.  Additional  tests  should  be 
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carried  out  for  the  MARINER  form  to  obtain  coefficients  Y6„  Md  Niu  t0 
prove  that  an  even  better  agreement  can  be  obtained  in  the  prediction  of 
tight  maneuvers  in  case  these  coefficients  are  included. 

Measurements  of  hydrodynamic  coefficients  should  be  obtained  for 
more  hull  forms,  for  which  reliable  full-scale  trials  are  available  so 
as  to  permit  comparison  with  computer  predictions  and  to  obtain  experience 
with  respect  to  prediction  accuracy. 

Corrections  of  hydrodynamic  coefficients  for  Reynolds  number  effect 
should  be  explored  further,  as  present  experience  is  insufficient  to 
permit  the  introduction  of  reliable  corrections. 
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Figure  1  -  Example  of  Measurements  of  Force  and  Moment  as  Function 

of  Drift  Angle 


rudder 

txtcutt 


Figure  2  -  Definition  of  Turning  Circle  Parameters 
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arm 


Figure  3  -  Diagrammatic  Definition  of  the  Zig-Zag  Maneuver 


Figure  4  ~  Results  from  Spiral  Maneuver  Showing  Slope  of  r-8  Curve 

for  Stable,  Marginally  Stable,  and  Unstable  Ships 
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Figure  5  -  Change  of 

Propeller  Revolutions 
as  Function  of  Speed- 
Loss  in  Maneuvers  for 
Different  Types  of 
Power  Plants 


forward  velocity  In  Knots 


Figure  6  -  speed  In  knots 

Velocity  Turn  Entry  Transient  for  35-Deg  Rudder  Computed  for  Different 
Types  of  Power  Plants 


Figure  7  -  sp,sa  tn  knots 

Time  History  of  Velocity  in  Zig-Zag  Maneuver  Computed  for  Different 
Types  of  Power  Plants 
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Figure  8  -  Results  from  Spiral  Maneuver  as  Influenced  by  the 

Assumption  of  Power  Plant 


Figure  9  -  Results  from  Spiral  Maneuver  Calculations  Presented 

as  Reciprocal  of  Turning  Radius  in  Steady  State  versus 
Rudder  Deflection 
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APPENDIX  A 


DATA  FORMS  FOR  PREPARATION  OF  INPUT  DATA 
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David  Taylor  Model  Basin 

Input  Data  -  Program  XPMC 

Page  -  1 

Hydromechanics  Laboratory 

Prediction  of  Maneuvers  of  Surface  Ships 

SPECIFICATION  DAW  -  SPECIFICATION  OF  MANEUVERS  TO  BE  PREDICTED 


CARD  SP  1 


LTCP>0,if  turning  circle  parameters  are  to  be  computed. 

LTC  specify  number  of  different  rudder  angles  for  which 
a  turning  circle  trajectory  should  be  predicted. 

LZZ  specify  number  of  zig-zag  maneuvers  to  be  computed. 

LSX  >0,if  spiral  maneuver  io  to  be  computed. 

Turning  circle  trajectory,  zig-zag,  and  spiral  maneuver 
will  be  plotted  if  GRAPH  >0. 

LTEST>0,if  new  set  of  input-data  ia  to  be  read  when 
this  computation  is  executed.  _ 


LTCP  * 
ISC  - 
LZZ  - 
ISM  - 
GRAPH* 
LT£ST= 


t 

1 

£ 

1 

0 


Specification  of  Turning  Circle  Parameters:  (If  LTCP>0  only)  CARD  SP  2 


DTC1  is  smallest  positive  rudder  angle  in  degrees  for  • 
which  parameters  should  be  computed. 

DTCD  is  the  difference  between  rudder  angles  for  which 
parameters  should  be  computed. 

DT02  is  maximum  rudder  angle  in  degrees  for  which 
parameters  should  be  computed. 


DTC1 

DTCD 

DTC2 


5.0 


5.0 


40.0 


Specification  for  Calc1  "ation  of  Turning  Circle  Trajectory: 
(Card  to  be  punchv  i  if  LTC  0  only) 

Calculation  of  the  turning  circle  trajectory  will 
be  carried  out  for  ITC  different  rudder  angles,  which 
should  be  stated  in  the  column  to  the  right. 

The  rudder  angles  should  be  g-.on  in  degrees,  positive 
for  port  or  left  rudder,  negative  for  starboard  or 
right  rudder. 


CARD  SF  3 


DTC(l) 

DTC(2) 

DTC(3) 

DTC(4) 


3  5.0 


J 


Specification  ror  Zig-Zag  Maneuver: 

(Card  to  be  punched  if  IZZ  0  only) 

The  zig-zag  maneuver  will  be  predicted  for  LZZ  dif¬ 
ferent  rudder  anglea  giver,  as  the  DZZ-values  in  the 
column  to  the  right. 


CARD  SP  a 


DZZ(l)  - 

tmmmmmn* 

20.0 

DZZ(2)  * 

DZZ(3)  - 

DZZC4)  * 

Specification  of  Spiral  Maneuver: 

(Card  to  be  punched  if  LSM  0  only) 

The  execution  of  the  spiral  maneuver  is  defined 
by  the  seven  parameters  in  the  column  to  the 
right  as  defined  in  the  figure  below. 


CARD  SP  5 


DSM(l) 

DSM(2) 

DSM(3) 

DSM(4) 

DSMD(l) 

DSMD(2) 

DSMD(  3  ) 


25.0 


10.0 


- 10.0 


* 25.0 


5.0 


1.0 


5.0 
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FORMAT' 9F8. 5)  FORMAT(9P8.5)  FORMAT* 9F8.5 )  FORMAT* 9?8. 5)  P0RMAT(9I8) 


David  Taylor  Modal  Basin  Input  Data  -  Program  XPKG 

Hydromechanics  Laboratory  Prediction  of  Maneuvers  of  Surface  Ships 


Page 


SHIP  DATA  I  -  SPECIFICATION  OF  SHIP  PARTICULARS,  ETC. 


? 


Title  of  Computations  The  title  will  be  printed  as  a  heading  CARD  SD  1 

on  all  result  sheets.  Mar.  72  characters.  F0RMAT(12A6) 


j  HYA  MAA/NEK.  PAS  DICTION  ASSUMING  CONSTANT  PQ*f£R  ^ 

Model  Identification; 

CARD  SD  2 

IMNO  is  the  model  number. 

IMNO  = 

6 295  I 

IPNO  is  the  propeller  number. 

IPNO  = 

6137  8 

Principle  Ship  Data  in  condition  tested  ; 

Length  between  perpendiculars. 

Used  for  nondimensionalizing  of  coefficients 

CARD  SD  3 

ft 

Airp  = 

528.01 

Length  on  waterline 

ft 

ALWL*= 

520.  76 

Beam  moulded 

ft 

BMLD*= 

7 6.02 

Draft  at  forward  perpendicular 

ft 

DFP*  = 

22.50 

Draft  at  aft  perpendicular 

ft 

DAP*  - 

25.70 

L.C.G.  measured  from  origin  of  axis  system 
positive  forward 

Radius  of  gyration 

measured  with  reference  to  origin  of  axis  syst 
Displacement 

ft 

CG*  » 

- 12.1 * 

ft 

em 

tons 

RAD*  = 

117.92 

DISP*= 

16005.0 

Number  of  propellers 

PROP  = 

1.0 

CARD  SD  4 

Propeller  diameter 

ft 

DIAM  = 

22.0 

Pitch  ratio  at  0.7  diameter 

PITCH* 

0.96 

Projected  propeller  blade  area  ratio 

AREA*= 

0.52 

Wake  cosff.  corresponding  to  condition  of  ship 

WAKE  = 

0. 160 

Thrust  deduction  coeff.  for  ship 

TDC  = 

0.136 

Approach  speed  (define  the  initial  condition) 

knots 

SPEED= 

15.0 

Rudder  rate  of  deflection 

deg/ sec 

RATE  = 

3.0 

Timelag  in  rudder  system 

sec 

TLAG  = 

0.0 

Type  of  propulsion  system 

TYPE  = 

- 1.0 

TfPE>  0.0,  if  constant  torque  (DIESEL)  LnMnM.1 

TYPE— 1.0,  if  constant  power  (TURBINE) 

TYPE<-5.0.  If  constant  propeller  revolutions. 


*  Data  marked  by  asterisk  are  not  used  in  the  calculations.  They  are  requested 
only  in  order  to  obtain  a  complete  description  of  the  ship  in  the  condition 
tested. 
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■;  ^  ;n>«w 


F0RMAT(9F8.2)  F0RMAT(9F8.2)  FORMAT(9IS) 


David  Taylor  Modal  Basin 

Input  Data  -  Program  XPMC 

Page  -  3 

Hydromechanics  Laboratory 

Prediction  of  Maneuvers  of  Surface  Ships 

SHIP  DATA  II  -  SHIP  EFFECTIVE  HORSEPOWER  DATA  and  OPEN  WATER  PROPELLER  CURVES 


Number  of  points  describing  ship  EHP-data: 

NEHP  $18 

Number  of  data  points  describing  propeller  curves: 
NPC  $18 


CARD  SD  5  m 

■*"'  "•  "■»U  H 

NEHPJ  9  S 
NPC  4  /7  ""  '"l  i 


4 


CARP  SD  6.1 


CARD  SD  6.2 


The  ship  EHP  should  be  defined 
for  a  speed  range  covering  the 
values  to  be  encountered  during 
the  maneuvers. 

Speed  values  in  knots  should 

VS(1)  - 

7.0 

CM 

• 

CO 

VS(10)  - 

/S(2)  - 

8.0 

fc 

O' 

Eh 

vs(ii)  - 

VSO)  > 

9.0 

O 

VS(12 )  * 

be  stated  in  the  columns  to  the 
right  in  order  of  increasing 

VS(4)  - 

10.0 

Pm 

VS(13)  - 

values. 

A  set  of  NEHF  speed  values 
should  be  given  (NEHP  is  de¬ 
fined  above). 

Only  one  punched  card  to  be 
used  if  NEHP  ^  9. 

VS(5)  - 

11.0 

VS(14)  - 

VS(6)  - 

12.0 

VS(15)  - 

VS(7)  - 

13.0 

VS(16)  - 

VS(8)  « 

tbo 

VS(17 )  - 

VS(9)  - 

15.0 

VS(18)  - 

CARD  SD  7. 

i 

CARD  SD  7.! 

Ship  EHP  data  corresponding  to 
the  speed  values  above  ehould 
be  stated  in  the  columns  to 
the  right. 

The  EHP  values  should  be  pre- 

EHP(1)» 

286.0 

CM 

• 

CO 

EHP(IO). 

EHP(2). 

4-14-.0 

o> 

& 

EHP(ll)* 

EHP(3)- 

582.0 

o 

EH?(12)» 

pared  using  a  l'oughnesa  allow¬ 
ance  corresponding  to  the  ship 

EHP(4)- 

789.0 

Pm 

EHP(13)« 

condition  for  which  the  maneu¬ 
vers  are  to  be  predicted. 

EHP(5)- 

1060.0 

EHP(14 )* 

EHP(6)- 

1391.0 

EHP(15)« 

EHP(7)» 

1815.0 

EHP(16)» 

EHP(8)» 

2318.0 

EHP(17)» 

EHP(9)» 

2919.0 

EHP(18). 

54 


.  jm*.  ■* 


David  Tavlor  Model  Basin  Input  Lata  -  Program  XTT.1C 

Hydromechanics  Laboratory  Prediction  of  Vr nenyf  r-  of  fu.  ■  ce  Ships 


Page  -  4 


SHIF  DATA  II  -  CONTINUED 

CANT  it  ft, 

,1 

CARD  SD  8. 

Open  Water  Propeller  Curves; 

ADV ( 1 ) a 

j  0.46  I 

/*’* 

’J" 

• 

»r 

ADV(10)n 

|  0.90  1 

Propeller  thruut  coefficient-, , 
and  torque  coefficients,!^, .should 

APV' 

!  050 

IT 

Th 

APVM.lV- 

j  0.95  j 

be  defined  for  advance  coeffi¬ 
cients  covering  the  range  encoun¬ 

ADVO)  = 

j  0.55 

1 

o 

ADV(1?)« 

i  | 

tered  during  maneuvers. 

Advance  coefficients  ADV*  V/(r»-Dl 
should  be  stated  in  order  of  in¬ 
creasing  values  in  the  columns  to 
the  right. 

ADV ( 4 ) » 

j  0.60 

ADVil’)«= 

i  & 

APV( ft )  ■ 

\  0.65 

! 

! 

APVf 14 ) - 

j  1 

ADV l ft)- 

]  0.70 

AI  V  lift)  ej 

'  ! 

A  set  of  NFC  data  should  be  given. 

( NPC  is  defined  on  Page-T).  If 

.'J'V  i  ’)- 

1  °S±  1 

l 

A '  V.  1  *  )  | 

j 

i 

KFC>9,  two  cards  should  be  unco. 

AW  8)- 

j  o.oo 

{ 

i 

i 

APVH7  '»| 

f  j 

ADV, '' ) - 

\_0.85 

i 

APV  { 1ft )  «j 

•ARP  CP 

i 

CARP  SP  0. 

FORMAT ( 9F8.S  )  FORMAT (9F5. 8) ”'5  FORMAT ( 9F8. 5) 


pi? 


F0RMAT(9F3.2)  FORMAT( 9F8.2) 


APPENDIX  B 

SAMPLE  OF  COMPUTER  OUTPUT 
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PAGE  1 


HY-A  MARINER.  PREDICTION  ASSUMING  CONSTANT  POWER. 


PREDICTION  OF  STANDARD  MANEUVERS 


MODEL  NUMBER  6295 

PROP.  NUMBER  613? 


TABLE  OF  CONTENTS 


PRINCIPLE  SHIP  DATA  {IN  CONDITION  TESTED)  .  .  PAGE  I 

CALCULATION  OF  COEFFICIENTS  XU.XUU.AND  XUUU  .  PAGE  2 

COEFFICIENTS  AND  STABILITY  ROOTS  .  PAGE  3 

TURNING  CIRCLE  PARAMETERS  .  .....  PAGE  4 

TURNING  CIRCLE  FOR  35.0  DEG.  RUDDER  ....  PAGE  5 

ZIG-ZAG  MANEUVER ( S )  .  ......  PAGE  7 

SPIRAL  MANEUVER  PAGE  9 
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PRINCIPLE  SHIP  DATA 


LENGTH  BETWEEN  PERPENDICULARS  .  =  528.01  FT 

LENGTH  ON  WATER  LINE . .  520.76  FT 

MOULDED  BEAM . .  76.02  FT 

DRAFT  AT  F.P . .  22.50  FT 

DRAFT  AT  A •  . . . .  =  25.70  FT 

L.C.G . ....  =  -12.14  FT 

RADIUS  OF  GYRATION . .  117.92  FT 

DISPLACEMENT  .  =  16005.  TONS 

NUMBER  OF  PROPELLERS . .  1.0 

PROPELLER  DIAMETER  .  =  22.00  FT 

PITCH  RATIO  AT  0.7  R . =  0.96 

PROJ.  AREA  /  DISC  AREA . .  0.52 

REVOLUTIONS . .  68.6  RPM 

APPROACH  SPEED  ........  .  *  15.00  KNOTS 

RUDDER  RATE  ..........  .  *  3.00  DEG/SEC 

TIME  LAG  OF  RUDDER  SYSTEM . .  0.  SEC 
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PACE  2 


HY-A  MARINER.  PREDICTION  ASSUMING  CONSTANT  POWER . 


CALCULATION  OF  THE  NON-DIM.  COEFFICIENTS  XU,  XUU,  AND  XUUU 
ON  THE  BASIS  OF  SHIP  EHP-D AT  A  AND  OPEN-WATER  PROPELLER  CURVES 


X-FORCE  =  PROPELLER  THRUST*! 1  — T  >  -  RESISTANCE 

=  XU*DELU  ♦  XUU*DELU* *2  ♦  XUUU*DELU*«3  (NON-DIM) 
WHERE  OELU  =  (U-UOl/U  =  NQN-DI M. CHANGE  IN  FORWARD  SPEED 


IS  COMPUTED  FOR 

APPROACH  SPEEO  UO  =  IS. 00  KNOTS 

WAKE  COEFFICIENT  W  =  0.160 

THRUST  DEDUCTION  COEFFICIENT  T  =  0.136 

TURBINE  SHIP  -  PROPELLER  TORQUE  VARY  PROPORTIONAL  TO 

REVOLUTIONS  IN  -1.000  POWER  DURING  MANEUVRE 


OPEN  WATER 

PROPELLER 

CHARACTER 

ADVANCE 

THRUST 

TORQUE 

COEFF. 

COEFF. 

COEFF. 

J 

KT 

KQ 

0.450 

0.291 

0.0446 

0.500 

0.270 

0.0422 

0.550 

0.250 

0.0398 

0.600 

0.229 

0.0372 

0.650 

0.209 

0.0345 

0.700 

0.187 

0.0317 

0.750 

0.165 

0.0288 

0.800 

0.143 

0.0257 

0.850 

0.  1  18 

0.0224 

0.900 

0.094 

0.0191 

0.9S0 

0.069 

0.0156 

EHP— I NPUT— DAT  A 

PROPELLER 

SPEED 

EHP 

RESIST. 

THRUST 

TORQUE 

*10-5 

*10-5 

*10-5 

KNOTS 

LB 

LB 

LB*FT 

15.0 

2919. 

0.6342 

0.7340 

3.0446 

7.0 

286. 

0.1332 

1.1010 

3.7643 

8.0 

414. 

0.1687 

1.0521 

3.6777 

9.0 

582. 

0.2107 

1.0064 

3.5919 

10.0 

789. 

0.2571 

0.9'*23 

3.5031 

It  .0 

1060. 

0.3140 

0.9  70 

3.4130 

12.0 

1391  . 

0.3778 

0.8  22 

3.3226 

13.0 

1815. 

0.4550 

0.8,0  74 

3.2303 

14.0 

2318. 

0.5396 

0.7?;.  16 

3.1370 

15.0 

2919. 

0.6342 

0.7  J41 

3.0447 

COEFFICIENTS  (NON-OIM) 
XU  =  -120.0E-5 

XUU  =  4S.0E-5 

XUUU  =  -10.3E-5 


TICS  ( INPUT  DATA) 


REVS. 

X  DELU 

X 

X*1 0*5 

*10-5 

*10*5 

FAIRED 

RPM 

LB  0-OIM 

O-DIM 

0-DIM 

o 

• 

CO 

55.5 

0.8181  -1.143 

2) 1.23 

211.23 

56.8 

0.7404  -0.375 

146. 36 

146.32 

58.  1 

0.6588  -0.667 

102.90 

103.02 

59.6 

0.5743  -0.500 

72.66 

72.51 

61.2 

0.4782  -0.364 

5C.00 

50.06 

62.8 

0.3758  -0.250 

33.02 

32.95 

64.6 

0. 2598  -0.154 

19,45 

19.54 

66.5 

0.1357  -0.071 

8.76 

8.77 

68.6 

0.0000  0. 

0.00 

-0.03 

59 
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HY-A  MARINER.  PREDICTION  ASSUMING  CONSTANT  PQP*ER» 
COEFFICIENTS  (INPUT  DATA) 


THE  FOLLOVING  COEFFICIENTS  OF  THE  EQUATIONS 
OF  MOTION  ARE  DESCRIBED  IN  0 TMB  REPORT  NO.XXXX 


X  -  E  0 

U  A  T  l  0  N 

Y  - 

■  E  0 

U  A  T  ION 

N  -  E  0 

U  A  T  I  0  N 

VAR¬ 

COEFF¬ 

VAR¬ 

COEFF¬ 

VAR¬ 

COEFF¬ 

IABLE 

ICIENT 

IABLE 

ICIENT 

IABLE 

ICIENT 

X  UDOT 

040 . 06— 5 

Y 

VDOT 

l S46.0E-5 

N 

VDOT 

-22.7E-5 

Y 

ROUT 

-8.6E-5 

N 

ROOT 

82  t 9E— 5 

X  U 

-120. OE— S 

Y 

V 

-1  16C.4E-S 

N 

V 

-26 3.SE-5 

X  uu 

45.0E-5 

Y 

VVV 

-8078. 2E-S 

N 

VVV 

1636.  IE— S 

X  OUU 

-10.3c— 5 

Y 

VRR 

-0.  E-S 

N 

VRR 

-0.  E— 5 

X  vv 

-898.BE-5 

Y 

VDD 

-3.8E-5 

N 

VDD 

12.5E-5 

X  RR 

18.0E-5 

Y 

R 

-499.0E-S 

N 

R 

— 166. OE— 5 

X  DD 

-94.8E-5 

Y 

RRR 

-0.  E-S 

N 

RRR 

-0.  E-5 

X  VK 

798 . OE— S 

Y 

R  VV 

15356. 0E-5 

N 

RVV 

-5433. 0E-5 

X  VO 

93.2E-5 

Y 

ROD 

-0.  E— 5 

N 

ROD 

-0.  E-5 

X  RD 

-0.  E-S 

Y 

D 

277.9E-5 

N 

D 

-138.8E-5 

Y 

000 

-90.0E-5 

N 

ODD 

45. OE-5 

Y 

OVV 

l  1 89.6E-5 

N 

OVV 

-489.0E-5 

Y 

ORM 

-0.  E— 5 

N 

DRR 

-0.  E-S 

Y 

OU 

-0.  E— 5 

N 

DU 

-0.  E-S 

Y 

VRD 

-0.  E— 5 

N 

VRD 

-0.  E-5 

X  0 

— 0 • OE— 5 

Y 

0 

-3.6E-5 

N 

0 

2.8E-5 

Y 

OU 

-0.  E-S 

N 

OU 

-0.  E-5 

OF  MASS 

« 

LB*SEC**2/FT 

O-DIM. WITH 

RH0*LPP»*3/2 

X  AND  Y  FORCES 

= 

LB 

- 

RH0»LPP**2*U**2/2 

N  MOMENT 

s 

LB*FT 

- 

RH0*LPP**3«0»*2/2 

U  AND  V 

s 

FT/SEC 

- 

U 

R 

ss 

RAO  I ANS/SEC 

- 

U/LPP 

0 

= 

RADIANS 

UDOT  AND  VOOT 

= 

FI /SEC/SEC 

- 

U**2/LPP 

RDUT 

RAO  I ANS/SEC /SEC 

- 

U**2/LPP**2 

STABILITY  ROOTS 


SIGMA  1 
SIGMA  2 
SIGMA  3 


-0. I 779E-00 
-0.2686E  01 
-0. 1429E-00 


STABILITY  CRITERION  n  0.61146-05 

SLOPE  OF  R—D  CURVE  =  -O.1037E-OO  <  DEG/SEC )/DEG 


60 
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MY-A  MARINER.  PREDICTION  ASSUMING  CONSTANT  POWER. 


TURN 

r  n  g 

CIRC 

l  e 

P  A 

RAM 

E  T  E 

R  S 

RUD 

ADVANCE 

TRANS¬ 

MAX 

TACT 

T  I  ME 

FOR 

MAX 

STEADY 

STEADY 

FINAL 

ANGLE 

(90  OEG ) 

FER 

ADVANCE 

DI  AM 

HDG 

CHANGE 

TRANS¬ 

-  TURN 

DRIFT 

SPEED 

(90  DEG) 

(  90  ) 

(  180  ) 

FER 

RAD 

ANGLE 

DEG 

FT 

FT 

FT 

FT 

SEC 

SEC 

FT 

FT 

DEG 

KNOTS 

5.0 

-5.0 

4859 

3830 

-3637 

2846 

4866 

3838 

76  15 

604  1 

278 

220 

530 

426 

-7623 

6051 

3771  . 
3001  . 

-3.9 

4.6 

13.77 

13.28 

10.0 

-10.0 

3175 

2838 

-2321 

2072 

3185 
284  7 

50  33 
4509 

134 

165 

364 
32  9 

-5042 

4519 

2508. 

2249. 

-5.6 

6.  0 

12.59 

12.28 

15.0 

-15.0 

2574 

2389 

-1852 

1  707 

2585 

240  1 

4091 

3796 

151 

140 

306 

285 

-4102 

3808 

2044. 

1898. 

-6.6 

6.9 

11.72 

11.50 

20.0 

-20.0 

2249 

2129 

-1582 

1496 

2262 
214  1 

3567 

3367 

133 

126 

275 

261 

-3580 

3379 

1785. 

1685. 

-7.4 

7.6 

11.02 

10.85 

25.0 

-25.0 

2049 

1  959 

-  1428 
1350 

2061 

1972 

3228 

3076 

123 

117 

25  7 
246 

-3241 

3089 

1617. 

*541. 

-8.  0 

8.  2 

10.41 

10.28 

30.0 

—  30.0 

1914 

1842 

-1309 

1246 

1927 

1856 

2990 

2866 

116 

1  1  1 

246 

236 

-3003 

2681 

1498. 

1437. 

-8.5 

8.6 

9.89 
9.  78 

35.0 

-35.0 

1818 

1760 

-1214 

1  167 

1833 

1775 

2817 

2714 

1 1  1 

107 

239 

231 

-2832 

2727 

1412. 

1360. 

-8.  8 

9.0 

9.  42 
9.  33 

40.0 

-40.0 

1751 

1703 

-1147 

1118 

1768 

1717 

269  l 

2598 

108 

105 

236 

228 

-2706 

2613 

1349. 

1302. 

-9.2 

9.3 

9.0) 

8.93 

bl 
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HY-A  MAR  I  NEW®  PREDICTION  ASSUMING  CONSTANT  POWER. 


T  U 

R  N  I  N  G 

C  I  R 

C  L  E 

FOR  35. 

0  D  E  G  • 

R  U  D  D  E 

R 

TIME 

RUDDER 

ADVANCE 

TRANS¬ 

SPEED 

HEADING 

ANGULAR 

DRIFT 

AFTER 

ANGLE 

FER 

ANGLE 

VELOCI TV 

ANGLE 

EXECUTE 

SEC 

DEG 

FT 

FT 

KNOTS 

DEG 

DEG/SEC 

DEG 

0. 

1.2 

0. 

0. 

15.00 

0. 

0. 

0. 

10. 0 

31.2 

252.7 

.  1 

14.91 

-2.0 

-0,462 

-1.6 

20.0 

35.0 

501.2 

-10.2 

14,59 

-9.8 

-0.972 

-5.2 

30.0 

35.0 

740.4 

-47.6 

14.13 

-20.3 

-1.066 

-7.6 

AO  »  0 

35.0 

962.6 

-118.7 

13.60 

-30.6 

-0.994 

-8.7 

5  0.0 

35.0 

1  162.6 

-220.5 

13.09 

-40.2 

-0.923 

-9.0 

60.0 

35.0 

1337.8 

-347. 7 

12.64 

-49.1 

-0.875 

-9.1 

70.0 

35.0 

1466.7 

-494.9 

1J.25 

-57.7 

-0.842 

-9.2 

80.0 

35.0 

1608,9 

-657.3 

11.90 

-65.9 

-0.816 

-9.2 

90.0 

35.0 

1704.2 

-830.8 

11.61 

-74.0 

-0.795 

-9.  1 

100. 0 

35.0 

1773.0 

-1011.5 

11.35 

-81.8 

-0.7*7 

-9.1 

110. 0 

35.0 

1815.8 

-1  196.0 

11.13 

-89.5 

-0.761 

-9.  1 

120.0 

35.0 

1833.3 

-1380.9 

10.93 

-97.  1 

-0.748 

-9.  1 

130.0 

35 . 0 

1826.7 

—  1 563. 5 

10.76 

-104.5 

-0.736 

-9.  1 

1A0.0 

35.0 

1797.2 

-1741.1 

10.60 

-111.8 

-0.726 

-9.0 

150.0 

35,0 

1746.2 

-191  l  ,Z 

10.47 

-119.0 

-0.717 

-9.0 

160.0 

35.0 

1675,1 

-2071.6 

10.35 

-126.1 

-0.709 

-9,0 

170.0 

35.0 

1585.8 

-2220.5 

13.25 

-133.1 

-0.701 

-9.0 

180.0 

35.0 

1479.9 

-2356. 0 

10.  16 

-140. 1 

-0.695 

-9.0 

190.0 

35.0 

1359.4 

-2476. 7 

10,07 

-147.0 

-0.690 

-9.0 

200.0 

35.0 

1226.3 

-2581. 2 

10.00 

-153.9 

-0.685 

-9.0 

210.0 

35.0 

1082.6 

-2668.5 

9.94 

-160.7 

-0.680 

-8.9 

220.0 

35.0 

930.6 

-2737. 7 

9.88 

-167.5 

-0.676 

-8.9 

230.0 

35.0 

772.3 

-2788.2 

9.33 

-174.3 

-0.673 

-0.9 

240.0 

35.0 

610.0 

-2819.7 

9.78 

-181.0 

-0.670 

-8.9 

250.0 

35,0 

445.8 

-2832. 0 

9.74 

-187.6 

-0.667 

-8.9 

260.0 

35.0 

281.9 

-2825. 1 

9.71 

-194.3 

-0.665 

-8.9 

270.0 

35.0 

120.5 

-2799. 4 

9.68 

-200.9 

-0.662 

-8.9 

280.0 

35  o  0 

-36.4 

-2755. 3 

9.65 

-207.6 

-0.661 

-8.9 

290.0 

35.0 

-186.8 

-2693.7 

9.62 

-214.1 

-0.659 

-8.9 

300,0 

35.0 

-328.7 

-2615.4 

9.60 

-220.7 

-0.657 

-8.9 

310.0 

3S.0 

-460.5 

-2521.6 

9.58 

-227.3 

-0.656 

-8.9 

320,0 

35.0 

-580.4 

-2413. S 

9.56 

-233.8 

-0.655 

-B.9 

330.0 

35.0 

-687.0 

-2292.7 

9.55 

-240.4 

-0.654 

-8.9 

340.0 

35.0 

-779.1 

-2160. 7 

9.53 

-246.9 

-0.653 

-8.9 

b2 

-  m  •ZSfiZZZjrj  -p*°~  ' 

-C  'J*  V"7* 

■':yz:-rz-^,'zrry, .. .. 

■ - nr 

ijwr-r"* 
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HV-A  MARINER.  PREDICTION  ASSUMING  CONSTANT  POWER. 


T  U 

R  N  I  N  G 

C  l  R 

C  L  F  F 

0  R  35. 

0  D  E  G  . 

R  U  0  D  E 

R 

<  CONTINUED) 

TIME 

RUDDER 

ADVANCE 

TRANS¬ 

SPEED 

HEADING 

ANGULAR 

DRIFT 

AFTER 

ANGLE 

FER 

ANGLE 

VELOCITY 

ANGLE 

£ “ECUTE 

oEC 

OEG 

FT 

F  T 

KNOTS 

DEG 

DEG/SEC 

DEG 

350.0 

3b. 0 

-855.4 

-2019. 3 

9.52 

-253.4 

-0.652 

-8.9 

360.0 

35.0 

-9  15.  1 

-  18  70.3 

9.51 

-260-0 

-0.651 

-8.9 

370.0 

35.0 

-957.5 

-1  715.  7 

9.50 

-266.5 

-0.650 

-8.9 

380.0 

35.0 

982.0 

-  155  r.  S 

9.49 

-273.0 

-0.650 

-8.9 

390.0 

35.0 

-988.5 

-  1  39  7. t 

9.48 

-279.5 

-0.649 

-8.9 

400.0 

35.0 

-976.9 

-  123h.  2 

9.47 

-285.9 

-0.649 

-8.9 

410„0 

35.0 

-947.3 

-1081.2 

9.47 

-292.4 

-0.648 

-8.9 

420.0 

35.0 

-900. 3 

-928.6 

9.46 

-296.9 

-0. 646 

-8.9 

4  30.0 

35.0 

-836.4 

-7  82.4 

9.46 

-305.4 

-0.648 

-8.9 

440  «  0 

35.0 

-756.4 

—  644  •  £ 

9.45 

-311.9 

-0.647 

-8.9 

450.0 

35.0 

-661.4 

-5 16. 4 

9.45 

-318.3 

-0.647 

-8.9 

460.0 

JS.O 

-552.6 

-400.0 

9.44 

-324.8 

-0.647 

-8.9 

470.0 

35.0 

-431.5 

-296.5 

9.44 

-331.3 

—0.647 

-8.9 

480.0 

35.0 

-299.5 

-207. 4 

9.44 

-337.7 

—0.6' 6 

-8.9 

490.0 

35.0 

-158.4 

-13  1.8 

9.44 

-344.2 

-0.646 

-8.9 

500.0 

35.0 

-9.9 

-76.5 

9.43 

-350.7 

-0.646 

-8.9 

510.0 

35.0 

144.1 

-36.  3 

9.43 

-357. I 

-0.646 

-8.9 

520.0 

35.0 

301.5 

-13.7 

9.43 

-363.6 

-0.646 

-8.9 

530.0 

35.0 

460.5 

—  8.9 

9.43 

-370.0 

-0.646 

-8.8 

540.0 

35.0 

o  1  9 . 0 

0 

9.43 

-376.5 

-0.646 

-8.8 

550.0 

35.0 

775.0 

-52.9 

9.43 

-382.9 

-0.645 

-8.8 

560.0 

35.0 

926.  S 

-101.  1 

9.42 

-389.4 

-0.645 

-8.8 

570.0 

35.0 

1071.6 

-166.0 

9.42 

-395.9 

-0.645 

-8.8 

580.0 

35.0 

1208.5 

-2 46.8 

9.42 

-402.3 

-0.645 

-8.8 

590.0 

35.0 

1335.4 

-342. 5 

9.42 

-408.8 

-0.645 

-8.8 

600.0 

35.0 

1450.8 

-451. e 

9.42 

-4  15.2 

-0.645 

-8.8 

610.0 

35.0 

1553.  1 

-573. 3 

9.42 

-42  1 . 7 

-0.645 

-8.8 

620.0 

35.0 

1641.1 

-705. 6 

9.42 

-4  28.1 

-0.645 

-8.8 

630.0 

35.0 

1713.7 

-847.0 

9.42 

-434.6 

-0.645 

-8,8 

640.0 

35.0 

1770.0 

-995.6 

9.42 

-44  1.0 

-0.645 

-8.8 

650.0 

35.0 

1809.2 

-1 149.5 

9.42 

-447.5 

-0.645 

-8.8 
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HY-A  MARINER.  PREDICTION  ASSUMING  CONSTANT  POWER. 


20.0  DEG.-  20.0  D  E  G  .  Z 


TIME 

RUDDER 

ADVANCE 

TRANS¬ 

AFTER 

ANGLE 

FER 

EXECUTE 

SEC 

DEG 

FT 

FT 

0. 

1  .2 

0. 

0. 

10.0 

-20.0 

252,8 

0.3 

20.0 

-20. C 

503.8 

10.  0 

30.0 

-20.0 

749.9 

42.  0 

40.0 

-5.0 

985.7 

103.  0 

50.0 

20.0 

1207.9 

191.7 

60.0 

20.0 

1420.7 

294.7 

70.0 

20.0 

1631.5 

399.  7 

80.0 

20.0 

1845.3 

497.  7 

90.0 

20.0 

2064.5 

579. 8 

100.0 

20.0 

2288.6 

638.  8 

110.0 

20.0 

2514,8 

670. 7 

120.0 

20.0 

2739.5 

674.4 

130.0 

20,0 

2959.3 

650.  8 

140.0 

5.0 

3171.3 

601.0 

150.0 

-20.0 

3373.9 

528.  0 

160.0 

-20.0 

3570.6 

442.  5 

170.0 

-20.0 

3766.9 

354.5 

180.0 

-20.0 

3966.9 

272.2 

190.0 

-20.0 

4172.3 

203.5 

200.0 

-20.0 

4362.7 

155.6 

210.0 

-20.0 

4595. 2 

132.3 

220.0 

-20.0 

4806.5 

135.2 

230.0 

-20.0 

5013.3 

163.  8 

240.0 

4.0 

5212.8 

217.  1 

2S0.0 

20.0 

5404.  1 

291.0 

260.0 

20.0 

5591.2 

375.4 

270.0 

20.0 

5778.5 

462.4 

280.0 

20.0 

5969.3 

545.9 

290.0 

20.0 

6165.4 

619.2 

300.0 

20.0 

6367.0 

676.3 

310.0 

20.0 

6572.6 

712.9 

320.0 

20.0 

6779.5 

726.8 

330.0 

20.0 

6985.0 

717.2 

340.  0 

20.0 

7186.3 

684.7 

G  -  Z  A  G  M  A 

N  E  U  V  E  R 

SPEED  HEADING 

ANGLE 

ANGULAR 

VELOCITY 

DRIFT 

ANGLE 

KNOTS  DEG 

DEG/SEC 

DEG 

15.00 

0. 

0. 

0. 

14.95 

1.9 

0.418 

1.5 

14.82 

8.  1 

0.74  j 

4.0 

14.60 

16.3 

0.863 

6,0 

14.32 

24.7 

0.738 

6. 7 

14.08 

29.  1 

0.  1  72 

4.8 

13.97 

28.7 

-0.185 

2.  1 

13.95 

25.4 

-0.430 

-0.6 

13.92 

20.0 

-0.626 

-3.  1 

13.63 

13.  1 

-0.739 

-5.  1 

13.66 

5.5 

-0.767 

-6.3 

13.44 

-2.2 

-0.754 

-7.0 

13.22 

-9.6 

-0.733 

-7.3 

13.01 

—  16.0 

-0.715 

-7.5 

12.83 

-23.6 

-0.600 

-7.2 

12.72 

-27.2 

— Oo 129 

-5.  1 

12.71 

-26.6 

0.J.92 

-2.3 

12.78 

-23.3 

0.423 

0.4 

12.83 

-18.0 

0.615 

2.9 

12.83 

-11.1 

0.732 

5,0 

12.74 

-3.6 

0.765 

6.4 

12.61 

4.0 

0.755 

7.  1 

12.46 

U.S 

0.735 

7.5 

12.31 

18.7 

0.717 

7.6 

12.20 

25.1 

0.507 

7.0 

12.  14 

27.7 

0.080 

4.7 

12.19 

27.  1 

-0. 161 

2.2 

12.28 

24.5 

-0.348 

-0.2 

12.38 

20.  1 

-0.S12 

-2.4 

12.42 

14.3 

-0,628 

-  4.4 

12. 41 

7.6 

-0.682 

-5.8 

12.34 

0.7 

-0.691 

-6.6 

12.25 

-6.  1 

-0.682 

-7.1 

12.  15 

-12.9 

-0.670 

-7.3 

12.05 

-19.5 

-  0.660 

-7.4 
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HV-A  MARINER.  PREDICTION  ASSUMING  CONSTANT  POWER. 


20.0  0  E  G  .  -  20.0  DEG.  ZIG-ZAG  MANEUVER 

( CONT 1 NUEO ) 


TIME 

RUDDER 

ADVANCE 

TRANS¬ 

AFTER 

ANGLE 

FER 

EXECUTE 

SEC 

DEG 

FT 

FT 

350.0 

-7.0 

7381.2 

629.  7 

360.0 

-20.0 

7569.5 

556.7 

370.0 

-20.0 

7755.  1 

475.6 

380.0 

-20.0 

7942.4 

394.0 

390.0 

-20.0 

3134.2 

318.  a 

400.0 

-20.0 

3331.8 

257.3 

410.0 

-20.0 

3534.2 

215.  2 

420.0 

-20.0 

8738.8 

196.  1 

430.0 

-20.0 

8942. 6 

201. 1 

440.0 

-20.0 

9142.5 

230.  2 

450.0 

4.0 

9336.0 

282.6 

460.0 

20.0 

9522. 1 

354.6 

470.0 

20.0 

9704.4 

437.  0 

480.0 

20.0 

9387.2 

522.5 

490.0 

20.0 

10073.4 

605.0 

500.0 

20.0 

10264.9 

678.6 

510.0 

20.0 

10462.2 

737.3 

520.0 

20.0 

10663.9 

776.7 

530.0 

20.0 

10867.6 

794.3 

540.0 

20.0 

11070.6 

789.2 

550.0 

20.0 

11270.4 

761.7 

SPEED 

HEADING 

ANGULAR 

DRIFT 

ANGLE 

VELOCITY 

ANGLE 

KNOTS 

DEG 

DEG /SEC 

DEG 

11.93 

-25.2 

-0.418 

-6.6 

11.96 

-26.9 

-0.000 

-4.1 

12.04 

—  25*5 

0.247 

-1.5 

12.  16 

-21.9 

0.448 

1.0 

12.24 

-  16.5 

0.612 

3.4 

12.27 

-9.8 

0.707 

5.2 

12.23 

-2.5 

0.733 

6*5 

12.14 

4.8 

0.725 

7.  1 

12.04 

12.0 

0.708 

7.5 

11.93 

19.0 

0.695 

7.6 

11.85 

25.2 

0.497 

7.0 

11.82 

27.8 

0.086 

4.7 

11.89 

27.3 

-0.147 

2.3 

12.01 

24.8 

-0.328 

-0.0 

12.  12 

20.6 

-0.487 

-2.2 

12.  19 

15.1 

-0.606 

-4.2 

12.20 

8.6 

-0.666 

-5.6 

12.  16 

1.9 

-0.680 

-6.5 

12.09 

-4.9 

-0.673 

-7.0 

12.00 

-11.6 

-0.662 

-7.2 

11.92 

-18.  1 

-0.653 

-7.4 

65 


PAGE  9 


HY-A  MARINER.  PREDICTION  ASSUMING  CONSTANT  POWER. 


SPIRAL  MANEUVER 


t 


RUDOER 

STEADY 

TIME  TO 

SPEED 

DRIFT 

TURNING 

ANGLE 

RATE  OF 

REACH 

IN 

ANGLE  IN 

RADIUS 

CHANGE  OF 

STEADY 

STEADY 

STEADY 

STEADY 

HEADING 

STATE 

STATE 

STATE 

STATE 

OEG 

1 JG/SEC 

SEC 

KNOTS 

DEG 

FT 

25.0 

-0.623 

605.0 

10.41 

-8.0 

1616.8 

20.0 

-0.S96 

437.0 

11.00 

-7.4 

1 785.6 

15.0 

-0.553 

482.0 

11.70 

-6.6 

2044.6 

10. 0 

-0.485 

537.0 

12.57 

-5.6 

2E0e.!3 

9.0 

-0.465 

372.0 

12.77 

-5.3 

2653.6 

8.0 

-0.444 

385.0 

12.99 

-5.0 

2831.3 

r.o 

-0.418 

400.0 

13.22 

-4.7 

3055.9 

6.0 

-0.388 

418.0 

13.47 

-4.3 

3353.  1 

5.0 

-0.352 

440.0 

13.74 

-3.9 

3773.0 

4.0 

-0.306 

467.0 

14.04 

-3.3 

4431.9 

3.0 

-0.245 

504.0 

14.38 

-2.6 

5684.3 

2.0 

-0. 151 

555.0 

14.75 

-1.6 

9423.4 

*.0 

0.004 

567.0 

14.99 

-0.0 

361035.4 

-0. 

0.  161 

568.0 

14.78 

1  .6 

8896. 1 

-1.0 

0.254 

529.0 

14.42 

2.6 

5494.5 

-2.0 

0.316 

484.0 

14.09 

3.3 

4318.0 

-3.0 

0.361 

452.0 

13.79 

3.8 

3690.5 

-4.0 

0.398 

428.0 

13.53 

4.3 

3287.5 

-■5.0 

0.428 

408.0 

13,23 

4.6 

3000.7 

—6 . 0 

0.454 

391.0 

13.06 

S.  0 

2783.2 

-7.0 

0.476 

377.0 

12.85 

5.  3 

2610.6 

-e.o 

0.495 

364.0 

12.65 

5.5 

2469.2 

-9.0 

0.513 

353. 0 

12.46 

5.8 

2350.6 

-10.0 

0.528 

343.0 

12.28 

6*0 

2249.2 

-15.0 

0.586 

480.0 

11.50 

6.9 

1898.3 

-20.0 

0.622 

439.0 

10.85 

7.6 

1685,4 

-25.0 

0.645 

405.0 

10.28 

8.2 

1540.8 

-20.0 

0.621 

425.0 

10.83 

7.6 

1685.5 

-15.0 

0.585 

466.0 

11.48 

6.9 

1898.5 

-10.0 

0.527 

514.0 

12.26 

6.0 

2249.6 

-9.0 

0.512 

348.0 

12.44 

5.8 

2351 . 0 

-8.0 

0.494 

359.0 

12.62 

5.5 

2469. 7 

-7.0 

0.475 

371.0 

12.82 

5.3 

2611. 

-6.0 

0.453 

364.0 

13.03 

5.0 

2783.8 

-5.0 

0.427 

399.0 

13.26 

4.6 

3001.6 

IN 
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HV-A  MARINER.  PREDICTION  ASSUMING  CONSTANT  POWER. 

SPIRAL  MANEUVER  ( CONTINUED  1 


RUDDER 

STEADY 

TIME  TO 

ANGLE 

RATE  OF 

REACH 

CHANGE  OF 

STEADY 

HEADING 

STATE 

DEG 

OEG/SEC 

SEC 

O 

• 

< 

1 

0.397 

417.0 

-3.0 

0.361 

438.0 

-2.0 

0.315 

465.0 

-1.0 

0.253 

502.0 

0. 

0.160 

554.0 

1.0 

0.006 

585.0 

2.0 

-0. 152 

567.0 

3.0 

-0.245 

531.0 

4.0 

-0.307 

486.0 

S.O 

-0.353 

454.0 

6.0 

-0.389 

429.0 

7.0 

-0.419 

409.0 

8.0 

-0.445 

393.0 

9.0 

-0.466 

378.0 

10.0 

-0.486 

365.0 

15.0 

-0.S55 

500.0 

20.0 

-0.597 

452.0 

25.0 

-0.623 

415.0 

SPEED 

DRIFT 

TURNING 

IN 

ANGLE  IN 

RADIUS  IN 

STEADY 

STEADY 

STEADY 

STATE 

STATE 

STATE 

KNOTS 

DEG 

FT 

13.50 

4.3 

3288.6 

1  3.77 

3.8 

3692. 1 

14.06 

3.  3 

4320.6 

14.39 

2.6 

5499.4 

14.75 

1.6 

8900.4 

14c  98 

-0.0 

250851.9 

14.77 

-1.6 

9419.1 

14.41 

-2.6 

5678.9 

14.07 

-3.3 

4429. 1 

13.77 

-3.9 

3771.2 

13.49 

-4.3 

3351.9 

13.24 

-4.7 

3055.0 

13.01 

-5.0 

2830.5 

12.80 

-5.3 

2653. 0 

12.59 

-5.6 

2508.0 

11.72 

-6.6 

2044.3 

11.02 

-7.4 

1 785.5 

10.41 

-8.0 

1616.8 
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INITIAL  CIECUTC  IK  IEC 


ZIC-ZAC  IW£WER 


STMIOMD 


•UDDER  AN«LC  U4  DC* 


l*0»7 


RATE  CP  CHANCE  OF  KAOING  VERSUS  RUDDER  ANCLE  FOR  PREDICTED  SPIRAL  MANEUVER 


APPENDIX  C 


FORTRAN  LISTING  OF  COMPUTER  PROGRAM 

The  computer  program  is  coded  in  the  FORTRAN  II  language  available 

for  the  IBM  7050  computer  at  IMB.  The  FORTRAN  listing  of  the  program  is 

included  on  the  following  pages.  Ih^  source  program  also  refers  to  the 

subroutines  AR  PLN1  and  AR  NXN1  for  the  least  squares  curve  fitting  used 

in  the  calculation  of  the  X  ,  X  ,  and  X  coefficients  and  AM  PLOT  for 

u  uu  uuu 

the  Charactron  Microfilm  Recorder. 

The  storage  required  by  the  program  can  be  greatly  reduced  in  the  case 

the  on-line  microfilm  plotting  is  left  out.  A  further  reduction  of  program 

length  and  storage  requirement  can  be  obtained  if  the  calculation  of  the 

coefficients  X  ,  X  ,  and  X  on  the  basis  of  EHP-ship  data  and  open- 
u  uu  uuu 

water  propeller  curves  is  carried  out  by  a  separate  program.  With  these 
reductions  of  the  program,  it  should  be  possible  to  run  the  program  on  any 
medium-size  computer. 

The  computation  time  for  a  prediction  of  the  "Standard  Maneuvers"  is 
approximately  6  min  on  the  IBM  7090 . 

Variables  in  the  program  have  as  far  as  possible  been  assigned  names 
that  correspond  to  the  established  nomenclature.  Tables  1-3  give  the 
relationship  between  the  hydrodynamics  derivatives  in  the  mathematical 
model,  Equations  (10),  arid  the  corresponding  identifiers  in  the  program. 


READ  5  •  QVDUT  •  QRD01  .  OV  ,  QV  VV  ,  GVR«  .  Q  VDD  ,  QM  ,  QRHR  .  OR  V  V.ORDB  .00*0000 


r  ui  x 

u  o  <  ♦ 

<  I  -  * 

o  k  a  * 
a  o  \  z 
a  o  —  < 
a  a  il  n. 
«  x  o 
a  <  • 
►-<*<■* 

<  jj  i  u 

•2*0 


0  3--  ~  I 

U  O  U.  -  u 

u.  a  *  a 

u  </)  a  rg  k 

0  0*1*  IX 

zoos*  •  <  ' 

<  X  ui  *  -  -  - 

>  ►■ujauvu  ( 

o  u  a.  o  *  \  *  < 

<  x  tn  a  n  w  «  i 

-  a  o  u  u  t 

•  •>  II  II  II  II 

~  -  n  n  «  m  ( 

**  U  O  O  O  U  C 


<  Z  uj 

I  >  o 

—  —  *  + 

■>  -»  a  n 

—  w  <  < 

‘  >  >  j  >  : 

J  QQUIQl 

<  <  o  <  ; 

•  II  n  n  II  : 

o  o  tf)  o  (/>  i 

!!<<<<( 
I  2  J  >  J  >  . 

Ml  IU  Q  "J  O  ( 

■  ->  o  <  o  «  J 


Ifl  —  J3 

<  CJ  < 

•  J)  I  J) 

>  o  a  ">  a 

:  o  x  -  a 

•  "  *  a  * 

)«/)(/>  U  /) 

;  <  <  -  < 

i  1/1  !•  |L  O 

■  a  o  *-  u 

I  X  *  z  • 

ii  «  -  m 

I  1/1  U  IIU 
:<  ii  i/i  ii 

i  z  <  J\ 

.  Q.  <  O  < 

:  a  i-  u  o 


i-  o  in 

•  •  <g  —  — 

n.  z  <  a  * 

c  a/  n  ~  x  u.  u 

a  >  u  a  a 
a  «  -*  *  z  in  z 

O  O  it/  *  *  • 

a.  •  a  —  — 

o  o  in  >-  ii  -»  ti 

•  -  ♦  h  « 

u  uj  ’**•  %n 

-j  a  a  o  >  o 

<  i  >  %  eu  ii  «o 

U  Ul  K  w  IIWI 

Wi  >  l« 


V)  •  «x 

-  J  a 

<  - 

o  a  a 


~  o  <v  a.  o 

O  •  —  •  ff>  c  «■» 

<  -  *o  "0  a  » 

fM  x  -  u.  in  u.  x  , 

-  -  O  01  o  " 


<«<«<«<< 

X  X  2  1  t  Z  7  1 

xxxaacz'xx 

OOOOCOOO' 

u.au.u.u.u.u.u.1 

-*  ♦  A  O  * 


J  l/l  3  !/»  «  «  ^  N  ( 

•>  7  >  B  J  X  X  . 

-  *  I  . 

4  a  —  x  —  -  wi  x  « 

*  <  i  s  Q  >*  — 

’*<*•<  >D  •  *  ; 

3  M  —  « 

’  O  Z  O  Z  X  >  -  ( 

J  Z  U  Z  3  •  •  -  V 

£ 

)  »  IA  >  (0  >>  |  C 

!  O  2  O  2  -  -  /\j  C 

/  J  a.  j  It  -  u.  X 

;  J  l  J  X  I  >•  >»  'I 

.  a  **  o  -  a  2  \  a 

IU,Oll,Q2«*i 


</)  uj  v  < 
a  a  i 

uj  u.  <  < 

>  a  uj 

3  H 

it-  z  in  < 

Z  O  < 

<  -  o  . 

*  »-  Z  t 

3  **  < 

>  Z  J  > 

u.  vi  a  t 

►-r  *—  UI  •— 

uu  u  • 

UJ  •  Z  J  < 

a  a  ~  a  < 

w  t/i  v)  c 

-j  J  c 

W  <  l/j 

O  >  ►*  I 

~  •  a  **  c 

a  u  z  . 

•-  j  ►-  -  j 

a  uj  z  j  c 

<  o  - 

a  •  i 

*“  uj 

<  j  s  a 

►-  UJ  M.  in 

<  Q  I*  J 

a  •  ‘n 

i  o 

*-  • 

D  O  u-  X 

a  <  j  j 

z  UJ  UJ  u 

-  x  o  a 


uoo  - 

n  n  *  oo  ii  o 

s  •  o  ^  * 

3  N  IJ  •  •  • 

5  n  I-  «  UI  . 

3  •  O  M  ^  ■ 

3  S-  »  •  w  •  < 

3  .T  -  •*  U  .*>  ' 

>  U  M  t.  <f  i 

3  -  *  a  < 

lao-'-o. 
)  U  •  U  •  N 

i  ►-  <0  *-  o  »g  > 


J3^>N.'0/I>0  OO 


•  j  axu^  —  x  • 

ot-  iiuauu*> 
-i*  iLzzaaav 
zw  z  *  *  z  z  2  • 

2k-  .  ,Xh 

.<  —II  II  «  -  «  O 

o  .j  a  ii  -3  ->  ii  ir  >  o 

z*  *ui  *  •-»-3>a 
a  u  o  a  •  «*  —  •  *x> 

-<ujz-in«  -  .  • 

♦  •  llJ 

aaai-i>-^  i  j 
■zavui/irohOijo 
■aj  *u/>ui<uud> 
■-<UZ---wwX>- 

>*«o  •»•*»»•• 

1O0UOCOQOOOO 

IUa/<U.liJUJ<lJUJU.u/Ut 

:aa:*axaaa£Z2a 


72 


O  UI 
Z  <  O 
a  <  ui 

w  a  o 


2  u.  a 

r  “j  n  . 

*“  u.  1 
O  1  Cl 

Z  Ul  - 

J  2  _!  I 

I  In  D  h  ■ 

-  n  J  O  u. 

X  4-  I  i  <  ' 

n  ^  i  <r  , 

<  c  n  a  ■ 

x  <  i 

X  i 

—  -  «r  ; 


tv  n  -  <  ~  cv  . 

—  m  -  a  -•  n  • 

<f  Cl  IV  “  Z  I  N  - 

ii  u  •-  ii  a  i 

tf.j2Xno.jZ: 

a  u.  a  o  •  a  u.  a  t 

2»atUl2"li 


.  *v  o  -«  a 

I  M  *.  \9 

-J  N  <  2 
♦  •  a  •» 

-  <">  z  r« 

v>  N  o  X 
<  (V  «V  -*  > 

>  a  it  —  i 

•  2  >'  \l  ~  9-  ! 

II  i- 

*  «  V1  y.  <  «  h  I 

•»  j  z  a  -•  2 

'  ~  y  —  —  , 

i  •  a  u.  ac  •  a:  < 
,  j.«  z  *-  a  <.  ai  a  ( 

i  o  -  x  o 

so  »n  >'■ 

i  m  rv  na  tv 


•*  X  «  • 
x  o  a 

N  «  Z  . 
•  •  x 
i  •  a 
i  x  a  uj 

•*  “  «-  i 

w  w  i  — 
♦“  *-  J 
<  <  2 
1  lul 
x  r  x  *j 
3  C  *  K 
a.  u.  ►-  < 
-  <V 

>0  O 
N  CD 

tv  *v 


*-  u.  2 

IX  1  3  Ui  L 

«  O  2.  O 

r  «  iu 

a  j  j  o  X 

I  •  •»  «t  w 

J  C  J  1 

<?  x  <  a  x 

X  •'*  X  -9  a 

-  «■  X  ~  2 

-  X  N  Q  I 

s  -  *  r  x 

i-  -  x  n  a 

u.  v  ->  <■  x 

*-  -«  x  « 

X  U.  N  —  - 

"  *-  —  X 

*2  Ifc  \  \ 

pj  tn  y»  i/i 

•  •  X  u.  2 


n  —  >  a  x  «t 

•  \  VI  <  I  X 

<v  cu  a  I  n  — 

• 

N  ■'  N  f  X  \ 

«  X  -  U  -v 

*  x  —  -  a»  U 

\  I/I  m 

•  •  H  \  (J1  \  V) 

\  N  H*  O 

*  *  tV  2  3  a/  X 

•  aza< 

*  •  k  a:  x 

U.  X  'V 

•  •  II  x  «S>  • 

I  o  ^ 

•  •  •  ift  •  CM  U» 


-»  a  i-  < 

3  iD  -t>  <  _) 

x  2  o: 

a  &  x  jj 

,  x  ~  j  t  a 

.  •*  3  r  3  ►« 

x  j  x  c  r 

-  c*  a  3  n 

r  \  n  «  ai  O 


D  ♦  0.  o 
O  2  * 
3  «  •  C 

/  a  \  <t 
~  2 

2  >-’»  Z  2 

*t  asa 
a  2  a  a 


N  ~  O  IV 

I  «  *  I 

—  u  s  ~ 

w  x  »W  — 

>1/1  *  > 

O  >  o  O 

<  >  (0  < 

•  Q  lO  ■ 

<  **  - 

■  ~  •  •  «- 

—  I/'  o  — 

3  <  T 

U  ■*'  >0  u 

—  a  ~ 


■  —  *  (V  I 

>0*0  I 

I  U  II)  ♦ 

•  •  <  w  I 

•  ~  W9  <  I 

i  m  a  v) 

i  a  a 

.  «  *  a  i 

•  —  l/)  • 

■  O  >  |T,  (/>  I 

i  \J  w  •>  <  | 

-  -  o  o  x  o  : 

,  u.  »  u  a  a  i 

•  N  ♦  *  CCN- 

;  2  u  m  ii  « 

•  -  it  o  -i  i 

ii  ^  ii  -»  o  < 

i  o  9  «  «  )• 

•  >  «-  *9  a 

■  o  x  ^  a  o 

i  y  ^  o  s  o 


o  < 

k  a 
ui  a 

u  9  o  a 
<  o  u.  x 
2  a  ui 
o  oz 

K  2  O  - 

a  <  n 

o  t-  it  *■> 

a  in  ui 
a  z  a  o 
a  o  >  o 
a  u  n 


•  u  < 

*-*  a  n  • 

U.  Z  ->  o 

a  •  «  y)  <o 

l/)  w  <  o 

*  •  noon 

u  i  » 

no  *  w  «  u 

—  \j  o  >  x  z 

m  *  <  u  o  ui  > 

>  in  o»—  ii  z  mu 

~un  >i  »>-oa 

ii  n  o  ui  —  j  x 

o^o>zaa  n 
uuoo>«>-- 


*-  fti 

,  —  *  tv 

h-  *  * 

u  .9  » 

>  •  >  in 

^  > 
l  fti  ,i  vA  • 

p  i  f  a  a 

o  «  «  a  <o 

>  -  •  • 

i  >  ►•  -rt  O' 

:  o  u  »•«  >  > 

i  <  -  i.  o  n 

s  a  p-  u  a 

.  o  k  ••  *  a 

p  u  z  to  n 

.  It  x  I  u  - 

i  in  n  «  n  n 

i  >  J'  i  ~  ~- 

•  y>  >  *  x 

.  a  ■  •:  «  a 

•  a  u  ■  o  a 


»n  u  u:  a 

3  O  -  V) 

x  a  w  v  i 

2X1 

a  UJ  »  “> 

•  z  -  — 

a  •  -  yi  n 

u  -*  •  >  > 

an^3  — 
a  n  n  x  ii 

*-  *  ii  « 

i-  o  r  -  n 

<v  w  “> 

x  oo  n  —  in 

n  t/i  o  > 

>  o  _j 
»  O  D  3  UJ 

>QX  XQ 


O  Q.  ~  ~ 

a  •*.  <v  —  a 

o  «  a  n 

_)  x  a  »»  3 

n  -»  _j  it  n  3  j 

II  II  <  o  3  3  3 

“j  M  o  X  X  X  X 


73 


PRINT  2.NPA&.TIT  4040  FOIMATC1H-.  7X , 59HC ALCUL AT ' ON  OF  THF  NON-DIW-  COEFFICIFNrS  XU.  XU 

PRlNTl/O.l^NO.lPNO  1U.  AND  XUUu/SX . 6 1 HCN  THE  OASIS  OF  SHIP  EHP-0ATA  AND  0P£N-hA T&W  PRO 

170  FORMAT ( 1 H— •  7X.63HP  RFOICTluN  OF  STANOArD  M  2P£LLtW  CURVF S/// U X . 4SHX-F0RCE  =  PROPELLER  THRUST»C|-T)  -  RESIbTAN 

1A  N  £  U  V  «.  R  ^////1  IX  I3HMC0EL  NUMBER  I  8/ 1  I  X 1  2HPRQP*  NUMBER  I  <*/////  3LF  l  )K  ,  5011-  XU«DELU  +  <U’-<«CjrLU»  *2  ♦  XLUU«OF  LU»  *  J  .  iVUN-0  I  M )  / 1  I  X  .  5 


in  — 

i  • 

iu  a  ® 

x  a  u. 

©  >  m 

-  I 

•  >  UJ 

-  X  » 

•  C4  * 

® 


I  >  CO 

>  u. 

•  >  ©  r 

i  i  - 

:  >  u 
x  >  v 

©  >  in 

•  xu 


X  •  ©  G 


lA  • 

i  «r  co 

u  x  u. 

x  z  © 

n  i 

•  z  uJ  Q 

-  X  Q 

•  ©  a  in 

CO 


I 

*  >■  UJ 


-  X  C  u'i 


i  ©  x  ♦ 

:  i  -  © 

■  «u  g  -  u. 

•23V 
i  ©  3  «n 


IT  -  - 


3  © 

3  I  Z 

■  Jj  X 

X  I  >  on 


Z  UJ  >  © 

X  >  u. 

©  a  © 

*  >  UJ 

-  x  a  • 

•  ©  3  i 

©  - 


•  >■  UJ  — 

-  x  o  v 

•  X  >  © 


x  > 

N  > 

-  o 


x  o  in  u 

H  I  UI 
*  Z  iu  in 


I  Q  © 

u  a  ul 

x  >  © 
n  i 
•  Z  Ui 
-  X  3 

x  •  ©  o 

N  ©  — 

n  u.  •  > 

v  « 


Ul  3  1 

x  g  n 


•  j 

a.  •  : 
a  o 
j  x 

•  a 

o 


a 

a  v- 
j  a 

V  Q 
3  > 


I  <  i 


3  O 

o  a 
•  a  x 

3  © 

X  © 

I  ©  — 

©  v 


U  U  MM 

UJ  ©  X  —  * 

©©MV* 

vs-  a 

3  3  ©  v  a 

u.  u.  z  ©  a  -j 

*  <  Z  0.  V 

©ll  -  <  J  M 

J  G  -  **  * 

<  O  M  * 

I  H  X  <  *  3 


M  Z 


v  a 
v  u 

V  • 


© 

Z  <7 

:  X  x 
muz 

:  ©  —  X  i 

x  v  m  : 


■VI  - 

—  a  « 
u  u  o 


j  ’j 
X  -  UJ 

X  ©  © 
3  <  V 
-1-0 
©  ©  u. 
X  X  O 
O  “  — 

—  V 
XXX 
M  <1 

-  V  - 
V  V  • 


© 


3  3  o  u 


o 
a  - 
o  © 


<3  - 

z  o 
—  > 

3  • 


©  o 
►*  x  l 
o  a 


o  v 
z  ~ 
<  o 


3  \  3  •  U 


M  J 
<  -  © 
J  X 


z  o 
a  3 

3  «n 


X 

• 

z 

© 

• 

Z 

a 

• 

Z 

JJ 

© 

• 

X 

M 

M  M 

© 

■T 

© 

3 

* 

■ 

3 

Q 

w 

X  M 

3 

— 

e 

3 

3 

3 

•? 

X 

> 

X 

3  a. 

o 

-  V 

\ 

V 

z 

* 

Z 

> 

© 

© 

3 

© 

V 

2 

i— 

z 

J 

2 

< 

z 

© 

Z 

© 

X 

© 

o 

X 

X 

V  M 

UJ 

© 

□ 

3 

> 

o 

1 

X 

M 

< 

a 

• 

a 

U 

+ 

a 

o 

a 

X 

© 

© 

o 

© 

cr- 

N 

© 

B 

o 

M 

M  • 

© 

z 

o 

ti 

1 

© 

© 

V 

o 

X 

a 

<7 

G 

2 

a 

o 

o 

2 

© 

z 

r> 

> 

u< 

U. 

O' 

3 

n 

> 

n 

u 

3 

V  * 

V 

< 

G 

3 

z 

1 

1 

3 

3 

o 

u 

> 

o 

3 

a 

< 

o 

X 

z 

V 

o 

o 

a 

3 

M  3 

o 

< 

> 

a 

3 

X 

o 

M 

3 

Z 

j 

< 

u 

a 

M 

© 

> 

© 

a 

3 

o 

o 

© 

3 

•  * 

© 

G 

J 

> 

a 

t. 

1‘ 

© 

o 

< 

1 

3 

3 

z 

a 

< 

> 

© 

1 

o 

< 

o 

© 

l 

o 

3 

< 

■  m 

X 

© 

< 

3 

> 

> 

Jl 

l» 

II 

UJ 

< 

r 

UJ 

M 

z 

II 

11 

3 

-J 

n 

c 

a 

X 

3  O 

z 

3 

a 

G 

u 

X 

3 

1 

z 

a 

3  • 

M 

X 

V 

z 

o 

It 

o 

> 

M 

3 

a 

z 

a 

o 

o 

z 

II 

3 

M 

o 

Z 

z 

z 

z 

> 

.1 

mm 

X 

a 

© 

3 

s 

© 

I 

j 

o 

z 

•  » 

M 

J 

3 

3 

J 

m) 

3 

2 

Z 

M 

mi 

H 

< 

> 

o 

1 

• 

n 

z 

> 

o  z 

1 

m 

z 

o 

m  n. 

V 

II 

< 

ii 

II 

C 

Q 

X 

7 

< 

z 

UJ 

< 

a 

Z 

• 

u. 

z 

01 

c 

© 

7. 

uj  a 

♦ 

3 

w 

a 

Z 

a 

Urn 

•  a 

V 

o 

G 

0) 

o 

©  © 

© 

© 

M 

U 

o 

G 

O 

o 

© 

© 

o 

a 

z 

a 

a 

— 

— 

M 

n 

© 

© 

3 

© 

* 

— 

— 

M 

n 

<  ©  o  s 

© 

-  M 

<7 

© 

o 

— 

M 

©  O 

V  z  z 

VXD 

o  - 

Jj  I  K 
UJ  ©  U 

a  •  3 
©MO 
•  u 
n  n  q 

t  u 

<  •  K 

B  II  © 

a  3 

o  o  ar 

u  3  r 

x  *■ 
z  <  i 

-  x  © 

©  M 
u»  *M  • 

O  •  X 

Z  Q  - 

<  UJ  - 

X  ©  V 

o  a  n 

•  ©  • 

X  © 

-Xu. 
o  ■>  • 

I  <  II 

7  a 
a  c 
Z  1  * 

a  x 

■I  <  < 

X  X 

3  <  <7 

V  -  M 

O  X  3  n 

3  -  Z  • 

I  -  ©  ®  . 

3  V  -  u 

-  Z  O  •  • 

3  -  li  . 

■I  ©  © 

U.  ( 

O  O  u  ►  1 

3  ©  3  X  < 

ui  3  iJ  <J  ■ 

G  3  X 

a  .J  M  ■ 

XX- 
O  <  • 

O  B  3  . 

X  Z  ■ 


>  > 
a  3 
<  iu 
>  Z 


o  o 
Z  z 
□  - 
3  X 
3 

a  g 

UJ 

3  X 
J  3 
UJ  B 

a  a 
a  a 
x 

a  x 


3  •  O 

XXX- 
R  X  -  U 
X  K  -  © 

&  -  v  uj 

<7  ©  o  3 


C  — 

a. 

>  3  < 

-  z  x  - 
~  -  z  * 

U  Z  O  Ul 

-  a  ©  > 


©  — 
UJ  UJ 
—  3 


O  O  X  X  o  o 

<7  —  ©  —  <t  M 

M  UJ  -  —  uj  M  W 

<7  C.  <  —  X  <t  <7 

>  »-  3 

a  v-  -  <  u  — 

>•  -  Z  X  5  -  z 

-  X  z  - 


u  u. 

3  UJ 

o  o 

<  u 

3  X 

<  © 
r  <m 


3  -  U 

©  —  a 
J  V  z 

J  Ul  • 
uj  3  -  - 

a  o  v  n 

□  3  V  1 
Z  3  O  • 
at-  k  * 


x  ©  • 

B  3  -  - 

r  n  < 

Z  X  *  -  • 

uj  *->  -  O 


©  z 

££  D 
3  - 

.  X  3 

■  M  J 

n  3 
>  *  > 

X  UJ  o 

-  a  © 

-  x  M 

-  <r  <r 
u  — 

<  *  u 
X  X  Z 
3  0- 
CM3 

l.  v  a 


•  u 


o  u. 
n  u 

•  7 
X  < 
lw  > 


OX  X  O  I 

I  ©  —  ^  * 

-  M  M  M  U.  1 

H  X 
<  04  >n  < 

a  -  -  z  x 

Z  V  S  -  3 

o  -  •  a  o 

u  <  u.  a  a. 

—  M 


©  *  - 
UJ  X  V 

a  -  v 
m  ©  a 

»  o  o 
a  >c  z  i 
x  ©  x  o 

UJ  - 
V 

©  X 
Q  ♦  X  - 
UJ  O  *  o 

uj  -  N  | 

a  •  v  o 
©  x  a 
ui 

a  x 


z  © 

UJ  ♦ 

3  3  0 

3  3  —  0 

U  UJ  *  3  © 

a  o  x  < 

ana 
r  -  a 


—  x  i 

•  a  > 

<  X  ©  3  < 


3  ©  ©  U  s 

a  >  x  *  < 

i!  ill  Cl  I  ! 

-  Z  -  J 

I  •  i 

a  © 
r  uj  i  < 
3  3  0©: 
XO  -  Jl 

©  z  * 

M  3  « 


<J  -  3 

3  « 

-  <  I 

Z  X  •  © 

-  Z  3  I 

X  CO  ©  o  ! 

a  u  -  —  ; 

—  M  f> 


3  a 
3  © 

X  • 


“>  w 


©  > 
>  — 
3  V 

a  3  — 

r  o  u. 

u  *  a 
z  3  i 
•  3  V 
-  3  — 
II  X  “7 
3  — 

-  a 
o  —  x 

n  II  uj 
<7  C  — 
*7  —  >1 
<  © 
I  o  u  > 

■  o  X  X 


©  — 
>  “> 
z  — 
i  a 
-  x 

U  UJ 
Q  • 

i  ~i 


a  3 


©  x 

3  © 


-  V 
3  U  © 
3  U  I 
3  U  ai 

x  a  x 
•  on 

3  X  • 
3  M  - 


X  M 


3  X 
X  3 


O' 


—  o  u.  o  o 

3  O  X  ©  I  3 

w  r>  •  UJ  <7  —  3 

X  <7  -  3  <7  »'  X 

3  3  Z  3  1 

|l  —  w  m  3  <  © 

©  Z  O  3  z  s  • 

>  -  O  Z  -  Z  X 

3  a  3  3  *  o  - 

x  a  x  o  a  l  - 


o  z 
o  <  * 
z  a  m 


3  0  3 

3  3 

u»  3  3 

ui  a  uj 

X  3  I 

3  ©  |  u. 

X  3  U. 

n  a  >-  3 

■«  o 

X  uj  o 

Z)  X 

V  <7 
X 

V  ?. 

—  o 

<  -  z 

>33  1 

3  <  3  3  X 

"02  « 

>  -  > 

>  3  u. 

*333 

3  a.  x 

X  z  ©  < 

•  —  «  x 

3  —  X  3  M 

C  ©  3 

uj  O  V 

Z  V  V 

J  ©  1  u  Z 

•  z 

3  3  H  |  3 

G  Z  3  X  - 

Z  < 

>33  3 


G  U  uJ  J 
>  IV 

—  3  V  -J 


x  c 
n  , 

•  z  ; 


Ul  —  u 
3  V 
Q  V  3 
<30 

-  z  o 

3  >  ■ 


.1  u-  > 

U  3  U 
a  Zj  © 

■U  <  I  3 
C  -  Ul  3 
U  to© 
©  a  i 


c 

3  c, 

z  > 


-  <  —  - 
a  a  a  a 
a  z  u.  a 


G  U.  X  I 

>  ©  « 

>  Ul  3  o  Z 

•  Z  Z 

3  J  J.' 

3  -  - 

Q  O  O  Z 

3  X  -  3 

x  ^  a 

o  X  3  z  z 

o  »  o 

.  <7  —  0  3  G 

-  z 

■  3  <  o  3  — 

z  a  z  o 

IL>  Z 

a  u.  a  -  < 


x  > 
n  > 
•  > 
-  3 
X  •  • 

—  X  > 

—  3  > 


-  M  *  © 


'A3  « 

I  O  Z 

©  o  ©  a 

iNO? 


74 


o 


u  o  o 

>  Z  U. 

D  -  - 
U'OX 
2  <  <0 
4  UJ  * 

2  X  -  - 


K  o 

Z  UJ  » 

O  UJ  l 

-  u  ax 

V.  —  J)  u  >• 

X  4  0 

X  di  4 

x  tr  u 

zx  *- 


CS  in  r 

UJ  u.  <0 

h-  a  ~  • 
u  jj  r  - 
-  a.  -o  ^ 
o  in  •  • 

UJ  £  —  — 

a  <  •  - 
a  r  —  * 


*  o 

Ul  4 
X  Ul 
<0  X 
-  X 
I  X  *n  O 
2  in  ui 
~  x  a  ~ 


—  *-  -  X 

r  x  ^  ~ 
N  m  x  < 
o  eg  a  >■ 


z  j  fa  D  ; 

U.  J  >  u 

o  a 

j  z  •  _»  i 

_»  <  —  j 

<  >  <  i 

u  a  a:  u  « 

-  »y 


u  —  <r  **  i  fj  o  : 


»  ^  *  n  u 

t  *  *  -  di  ' 

.  O  X  X 

b  UJ  &  -O  d)  l 

(  on, 

'  >  -  «  _l  < 

»  ►-  < 

3  o  -  o  «  : 

:  cm  x  j  t.  . 

•  u  o  a  i  - 

-  w*)  <  2  2 

-2  Xu.-, 

too 

L  X  i 

■  J  o  ; 

3  x  o  u.  <  a  « 

L  O  -  O  U. 


Z  UJ  UJ 
-  >  > 

X  C  uJ 

o  o  o 

JUOZ 
J  - 
O  J  J  *- 
Z  _J  J  2 

<  <  <  a 

uuu 


♦  a  o  o  i  i  : 

i  o  2  o  n  x  i  • 

<  •  o  o  -  -  J 

CL  !\J  fU  IV  —  —  * 

2  J-  >-  I 

II  l-  u  *-  <  <  • 

-  o  z  z  z  x  x  : 

n  <  -  -  -  r  x  : 

-  cl  c  a  a  o  o  < 

x  2  a  a  a  u.  u.  \ 


►“  in  ~  N 

4  x  n  n 

|  N  «  «  111 


-  j)  eg  -  .*1  * 

x  *-  i 

V)  O  u  <  N  I 

-  o  u.  •-  z  x  ->  ; 

II  It  O 

c  u.  ii  o  «  o  n  : 

n  o  o  o  a  u.  -  t 


X  >•  x‘ 

o  o  u 
uj  <  z  a  u 

uj  x  uj  in 

>-  *-  X  u. 

in  u 
w 

o  - 

X  X  Q  O  »- 

4  X  «  u. 

♦  I 


I  X  O 

.  iu 

d.  o  — 

O  k-  s 

-  It  N 

:  —  r  x 

t  O  CD  X 

i  U  <A 

OX  U 

•r  j 

0*  0  2 
O  S  UJ  X 

)  —  o  o 

i  ‘-j  uj 

:  j  a  x 

•  O  Jl  C  J 

Z  *-  Uj 


4  • 

«-  N 

~  o  - 


2  Z  Uj 

i  —  a  -i 

i  1  3  O  U 


1  4  £  - 

jj<  oin 

>  <V  n  *  if* 


IlOO  •  - 

*  *M  Ml 
-  o  M  a.  - 
U  ,1  »  -  'Jl  f\|  N  1 
►*  b-  X  I 

oo»-4  i  o  h  - 
--t-ZX^XZ' 
li  ll  -  X  I'  -  — 

1/iu.oacou.X' 

xooau.*«-a.' 


u.  •  J 

C3  O  • 

o  - 

Z  N.  || 

o  —  -  — 


<  o  o  -* 

J  *-  N  u 

O  J  N  H 

U  —  n  -  o  • 


>  r  j  v*»  2  Z  : 

||  <  m  M  r 

.mainn 
'  x  z  a  a  c 


x  ; 
;  o  i  n  a  • 

-  N  < 

■  D  -  -  J 

*-  ►-  o  : 

:  x  <  <  o  » 

J  X  7  < 

:  •  r  x  <  i 

I  O  O  I 

,  o  u.  u. 


<  UJ  O  f\J  * 

ll  O  O  —  —  O 

iw  z  w  *  i-  a  o 

,  >  —  t  n  x  fd 

X  *-  n  “>  >.  ii  i> 

i  —  z  -  i«  -  -  *y 

■  u.  o  >•  n  u.  v  >• 

:  -  u  a  —  -  <x  a 


75 


PtflNT  *U?C*  IF((D-CjF  )  >  bZBG  ,‘,240 .6240 

JJ^JI-UVb  b24'J  D  =  DF 

IF  (  (OF  -ObVC  J  J  J  >  •b(»:>>  >.*<  40«^Oo*^.»OO  *TC  ^2 


Q  • 

a.  cs 

>  • 


o  * 

a  a 

a  • 

*  > 
>  * 

*  o 

>  a 

•  > 

>  a 


a  - 

•  3 

Z  J 


a  a 

—  • 

+  a 


—  a  > 


Oil 
I  X  •  ' 

cr  t 

■  a  a  < 

>  a  >  ■ 


i  O  *  : 
I  X  > 

!  X  *  , 
i  ♦  >  : 

x  >  ■ 

1  3  >  • 

i  •  >-  < 


u  u  i 
V  u  l 

<  <  • 


2  C 
+  u 
>  c 


■  K  X 

;  j  - 

l  d)  o 


o  — 
(  3l  - 

w  n 

a.  a 
.  »-  o  — 


—  3  X  t>  > 

l  I —  c  > 

-  a  —  >  — 

1  l  O  \  IM 
«-  t.  e  »  . 


O  n 

N»  N 

CM  N 
J>  Si 


o  o 
N  flO 
N  CM 
01  lT 


-»  O  >  —  < 

a  a  ■ 

•  0  ♦  —  O 

«  O  «  •  ( 

avast- 

:  o  *  «  o 

i  •  a  cn  •  ! 

O  O  O  O  i 

—  <M  —  < 


«-  <  ■ 
i  a  a  : 

)  Q  > 


3  3 
O  O 
3  > 


'  <  *  • 

I  N  H  - 
I  X  3  3 

-  o  o 

3  3  * 
I3  +  + 
*11 
■  a.  —  — 
-go 

•  3  3  > 
J  i.  « 

M  a  a 

I  O  -  - 
l  U  C  Q 
;  «  3  > 


•»  c 
a  a 


O'  o 
O  CM 
IP  u> 

•  in 


X  Ci 

-  a 

Q  «-* 


«  N  -  jJ  3  3 

t-  O  3  U  O  > 

3  0  3  »  •  3  J 

0  3  n  a  >  o  ~ 

r  i  >  t.  - - 3 

■*  X  —  0.  O  •  •  ♦ 


a.  —  j.  U.  v 


O  'V  - 
►-  O  — 

a  01  n 

•  •  a 

—  O  — 
>  N  vt 

-  n  * 


-  _  >  z  •  ^ 

■  *-  •  <  o  o 


-03 
3  3  S 

a  -  a 

it  I*  — 

x  3  a 
-  j  » 

O  oJ  II 

a  o  > 


<  ♦  J  J  -  -  3 

m  ►*  j  x  ^  3-1 

I  O  3  Q  >  X  ii 


I  u. 


<*•• 

i  iT  n  •  <n  nwuoo 

ii  -  a  •  x  -  i  i/i  x  >■  »- 

ii  ii  ii  a  u.  ►*  a  i  -  i  —  ii  -  m  n  ii  <i  » 

,  ii>x2ani»xi-'-  i*'iixx<ooo 
-JJ--ZO-  «  'Ml  u  a - ?  *  ’ 

i/lUuiOO  —  <  C  3  u.  U.  0)  t 

aoovxu.ao-3-aa  : 

o  o 


-  a 
o  - 

1X0 

I  I|  X  <f 


*■  -  a  -v  < 
a  3  «  -  - 


<  O  tf*  o 

—  X  'I  > 

■  3  -  U  — 


a-zxa-zza* 


>  o  • 

"  >■  «  *- 
i*  i'  n  2 
a  -  u  3 
v  a  i-  a 
z  >■  a  « 


3  3 
CM  3 
-O  »0 
n  o*> 


in  *o 
>o  *c 
n  01 


3  o 

<y  o 

s>  <n 


UJ  oj  Z  - 
JO<  • 

ooxo 

Z  3  U  - 
<  a  • 
i l  - 


a  o  a; 
oj  a  *- 
>  x  -n 
3 
O  O' 

z  a  X 
-  <  <r 
o  »-  in 

<  m  —  i 
uj  X  | 
I  n 
I  >0 


uj  x  a  ■ 

O  CD  o  ' 

z  o  a  i 


-  X  1 
O  N  . 
C  N  ■ 

J  X 

-  'M 
D  X 


a  - 
-  e 
o  a 

a  x  o 


3  O 
-J 

-  2 

■  \J  —  1 

I  O 


o  no 

o  -  ^ 

^  N  PI 

**J  cv  M 


3  -  a',  <  O  ^ 

■  .  a  j 

’5 

i  5  < 

u. 

•  a 

-  CM 

o  n 

K  ffl 
*M  CM 
CM  CM 


^  r»  a  o 

X  •  3  T 

*n  u 

u.  u.  3 


X  J 

kU  >C  < 

o  •  u 

*n  <r 


Si  3 
3  - 
3  k- 
3  3 

z  a 

-  z 

-  i  - 

z  -  o 

O  X  ll 

u  m  o 


a  a.  n 

-  c  -  \ 

i  o  >  •  CM 

--•O'- 

x  -  a  a.  in 

—  in  —  *  • 

u.  *  c  —  in 


<  « 
3  3- 
k  j  n 


n  o  o  • 

u  3  0 

a  uj  i.  -  ii 

n  -  ^  a  ii  a 

z  i*  <  —  u  ii 

—  «  a  a  n  in 

a  t  u  3  l> 


z  -  n 

*  ~ 

3  -  ^ 
3  J  | 
N.  UJ  — 

in  3  a 

ii  w 

o  -  u. 

Q  —  — 


•  O  - 

3  -  n 

•  u.  — 

—  •  c 

—  —  n 

o  S'  in 

cm  u.  i 

—  rJ 

O  W  f 

t-  — 

I"  <  X 

z  a  - 


01  •  oi  u. 

—  j  3  a 

p  -  <  —  • 

3  ♦  a  o  o  < 

i  o  z  —  ^  i 

<  •  3  <t  ■ 

^  -  a  ?m  i  -  , 
♦  —  z  - 
t-  <t  x  it  h.  a  h-  i 

z  ->  W  o  z  -  z  i 


I  cm  O  < 
K.O( 
*»  : 


C  X 
00  *M 
PJ  0) 

-  0*  * 

CM  *  O  * 

X  O  ®  t 

-  —  <M  i 

i  IM  .1  - 

a  -  n  -  . 

-  <  *  3  > 


CM  *  /I  I  ( 

x  »-  <->  a  « 

i  3  -  a  i 


<  —  —  —  —  • 


~  -  —  *  -  ' 
zzxx*<xn< 


1  -  H  —  -  3 


76 


«we»a>ip^ 


■»3P= 


DELSP)  S72I. 5721 .5722 


Ul  N  *  - 


U.  1  I  «  W 


10  I 
CD  O  O 

<  a  q 
n  d  a 
*•  <  < 
o  «  w 

o 

3ILIL 

<  -  - 


*  M 

UJ  • 

z  o 


x  a 

10  —  Ul  N 
K  UL  D  w 
Ik  CO  z 

I  CD  -•  □ 
*  <  K 

W  W  z 

11  U.  0  Q 
W  N  U  O 


i/>  N 
U  O 
K  O' 

J  ~ 


Zd 

UJ  Ik 


a  * 

UJ  to 
UJ  UJ 

a  h  <u 

Ifl  JOO 

*  ♦  o 

*  fO 

O  H. 
O  Q  »-  Z 

<  <  fr- 

iD  uj  o  a 
a  a  o  a 


w  o  ~ 
a  O'  o 
•  x 
no  a 
a  o  x 
*■>  —  I 
>  ♦ 
Q  *  • 
Z  O  « 
uj  e  . 


_  -  m 
~  a  a 
Iff  <  8  *) 
O  X 
4IUUIO 
*-  a.  o  j  z 
»■<  j««i 
<  cc  u. 

X  O  J  J 
E  -  J  O  J 
Olt<Z< 
u.  «  u  w  u 


-NiflOO 

N  N  w  n  * 

N  S  N  f>  N 

m  ui  ui  m  in 


o  o  o  o  o 

«  wn  in  ^ 

o  o  o  o  o 

O  O'  O  O'  o 


Q 

Z 

UJ 


77 


■wr^ 


REFERENCES 


1.  '‘Nomenclature  for  Treating  the  Motion  of  a  Submerged  Body  Through 
a  Fluid,"  The  Society  of  Naval  Architects  and  Marine  Engineers,  Technical 
and  Research  Bulletin  No.  1-5  (Apr  1950) 

2.  Imlay,  F.B.,  "A  Nomenclature  for  Stability  and  Control,"  David 
Taylor  Model  Basin  Report  1319  (May  1959) . 

3.  Gertler,  M.  and  Gover,  S.C.,  "Handling  Quality  Criteria  for  Sur¬ 
face  Ships,"  David  Taylor  Model  Basin  Report  1461  (Oct  1960). 

4.  Abkowitz,  M.A.,  "Lectures  on  Ship  Hydrodynamics  -  Steering  and 
Manoeuvrability,"  Hydro-  and  Aerodynamics  Laboratory,  Lyngby,  Denmark, 
Report  Hy-5  (May  1964) . 

5.  Chislett,  M.S.  and  Strom-Tejsen,  J.,  "Planar  Motion  Mechanism 
Tests  and  Full-Scale  Steering  and  Manoeuvring  Predictions  for  a  MARINER 
Class  Vessel,"  Hydro-  and  Aerodynamics  Laboratory,  Lyngby,  Denmark, 

Report  Hy-6  (Apr  1965) .  (Also  published  in  International  Shipbuilding 
Progress,  Vol  12,  No.  129  (May  1965)). 

6.  Strom-Tejsen,  J.,  and  Chislett,  M.S.:  "A  Model  Testing  Technique 
and  Method  of  Analysis  for  the  Prediction  of  Steering  and  Manoeuvring 
Qualities  of  Surface  Vessels,"  Hydro-  and  Aerodynamics  Laboratory,  Lyngby, 
Denmark,  Report  Hy-7  (under  preparation) 


78 


^HrSSa STfi.  ■  jr" 


;.ts 


UNCLASSIFIED 
Security  Classification 


DOCUMENT  CONTROL  DATA  -  R&D 

(Security  classification  of  titto,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  is  classified) 

1  ORIGINATIN  G  ACTIVITY  (Corporate  author) 

David  Taylor  Model  Basin 

2a  REPORT  SECURITY  C  L  ASSI  FIC  A  TION 

Unclassified 

Department  of  the  Na\*y 

Washington,  D.  C. 

2b  GROLP 

3  REPORT  TITLE 

A  DIGITAL  COMPUTER  TECHNIQUE  FOR  PREDICTION  OF 
SHIPS 

STANDARD  MANEUVERS  OF  SURFACE 

4  DESCRIPTIVE  NOTES  (Type  of  report  and  inclusive  dates) 

5  AUTHORfS.)  (Last  name,  first  name,  initial) 


Strom-Tejsen,  Jorgen 
_________ 

December  1965 


7a  TOTAL  NO  OF  PAGES  7b.  NO  OF  REFS 

87  6 


8a-  CONTRACT  OR  GRANT  NO. 

9a.  ORIGINATOR'S  REPORT  NUMOERfS.) 

b.  PROJECT  NO 

2130 

c 

9  b.  OTHER  REPORT  NO  (S)  (A  ny  other  numbers  that  may  be  assigned 

this  report,) 

10  AVAIL  A8IHTY/L1MITATION  NOTICES 

Distribution  of  this  document  is  unlimited 


11  SUPPL  EMENTARY  NOTES 


13  ABSTRACT 


12.  SPONSORING  MILITARY  ACTIVITY 

Department  of  the  Navy 
Bureau  of  Ships 
Washington,  D,  C. 


This  report  presents  a  computer  program  for  the  solution  of  a  mathematical  model 
representing  the  motion  of  a  surface  ship,  giving  predictions  of  steering  and 
maneuvering  qualities.  The  nonlinear  mathematical  model  based  on  a  third-orde. 
Taylor  expansion  of  forces  and  moments  in  the  equations  of  motion  is  reviewed . 
The  hydrodynamic  derivatives  representing  the  input  to  the  program  can  be 
obtained  from  present  captive  model  testing  techniques.  Any  motion  of  a  sur¬ 
face  ship  including  tight  maneuvers  and  loop  phenomenon  recognized  in  the  spiral 
maneuver  for  a  directional  unstable  ship  should  be  predictable  with  accuracy. 

The  computer  program  which  gives  predictions  for  the  "Standard  Maneuvers,"  turn¬ 
ing  circles,  zig-zag,  and  spiral  maneuver,  is  described,  and  results  of  sample 
calculations  are  included.  An  instruction  for  preparation  of  input  data  for 
the  program,  samples  of  the  computer  results,  and  the  FORTRAN  listing  of  the 
computer  program  are  also  given. 


FORM 

1  JAN  G4 


1473 


UNCLASSIFIED _ 

Security  Classification 


IKCHW 


ROUE 


Hydrodynamics 

Maneuvering 

Steering 

Computer  Program 


INSTRUCTIONS 


ORIGINATING  ACTIVITY:  Enter  the  name  and  address 
of  the  contractor,  subcontractor,  grantee,  Department  of  De¬ 
fense  activity  or  other  organization  ( corporate  author)  issuing 
the  report. 

2a.  REPORT  SECUMTY  CLASSIFICATION:  Enter  the  over¬ 
all  security  classification  of  the  report.  Indicate  whether 
“Restricted  Data”  is  included.  Marking  is  to  be  in  accord¬ 
ance  with  appropriate  security  regulations. 

2*.  GROUP;  Automatic  downgrading  is  specified  in  DoD  Di¬ 
rective  5200.10  and  Armed  Forces  Industrial  Manual.  Enter 
the  group  number.  Also,  when  applicable,  show  that  optional 
markings  have  been  used  for  Group  3  and  Group  4  as  author¬ 
ized. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letters.  Titles  in  all  cases  should  be  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  show  title  classification  in  all  capitals  irt  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.  g.,  interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  Enter  the  name(s)  of  authors)  as  shown  on 
or  in  the  report.  Entei  last  name,  first  name,  middle  initial. 

If  military,  show  rank  and  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  date  of  the  report  as  day, 
month,  year,  or  month,  year.  If  mire  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  i.e. ,  enter  the 
number  of  pages  containing  information. 

7b.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  report. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

8ft.  8c,  &  8 d.  PROJECT  NUMBER;  Enter  the  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a.  ORIGINATOR’S  REPORT  NUMBER(S):  Er.ter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  umque  to  this  report. 

9ft.  OTHER  REPORT  NUMBER(S):  If  tho  report  has  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  the  sponsor),  also  enter  this  number(s). 

10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 


imposed  by  security  classification,  using  standard  statemer. 
such  as: 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC” 

(2)  “Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized.  *’ 

(3)  "U.  S.  Government  agencies  mav  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 


(4)  "U.  S.  military  agencies  may  obtain  copies  of  th 

report  directly  from  DDC  Other  qualified  users 
shall  request  through 


(5)  "All  distribution  of  this  report  is  controlled.  Qual¬ 
ified  DDC  users  shall  request  through 


If  the  report  has  been  furnished  to  the  Office  of  Techmca 
Services,  Department  of  Commerce,  for  sale  to  the  public,  ind 
cate  this  fact  and  enter  the  price,  if  known. 

11.  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  na-,e  of 
the  departmental  project  office  or  laboratory  sponsoring  (pay¬ 
ing  forj  the  research  and  development.  Include  addr  -ss. 

13.  ABSTRACT:  Eater  an  abstract  giving  a  biiei  and  factua 
summary  of  the  document  indicative  of  the  report,  even  thoug'- 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re¬ 
port,  If  additional  space  is  required,  a  continuation  sheet  sh 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  rep< 
be  unclassified.  Each  paragraph  of  the  abstract  shall  end  wit 
an  indication  of  the  military  security  classification  of  the  in¬ 
formation  in  the  paragraph,  represented  as  (TS).  (S).  (C).  or  (l 

There  is  no  limitation  on  the  length  of  the  abstract.  How 
ever,  the  suggested  length  is  from  150  to  225  words. 

14.  KEY  WORDS:  Key  words  are  technically  meaningful  term 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Key  words  must  be 
selected  so  that  no  security  classification  is  required.  Identi 
fiers,  such  as  equipment  model  designation,  trade  name,  militi 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rules,  and  weights  is  optional 


UNCLASSIFIED _ 

Security  Classification 


or  -  j&sKmaszr.  - 


